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Abstract. The spatial gradient of solutions to nonlinear degenerate parabolic 
equations can be pointwise estimated by the caloric Riesz potential of the right 
hand side datum, exactly as in the case of the heat equation. Heat kernels type 
estimates persist in the nonlinear case. 
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1. Results 

In this paper we are going to consider nonhnear, possibly degenerate parabohc 
equations whose model is given by 



(1.1) ut - div {\Du\P-^ Du) ^ 



where fi denotes in the most general case a Borel measure with finite total mass. 
Although the kind of problems considered are nonlinear our goal is to provide a 
suitable linear potential theory aimed at describing, in a sharp way, the regularity 
properties of the gradient Du in terms of those of fi. More precisely, our description 
shows that sharp gradient pointwise estimates can be given in terms of classical 
Riesz caloric potentials of the right hand side fi. We will see that, surprisingly 
enough, bounds similar to those that hold for the heat equation 

(1.2) Mt - Alt = /X 

actually hold for solutions to ()l.ip and, more in general, for solutions to quasilinear 
equations of the type 

(1.3) ut - div a{x,t, Du) = fi in = fl x {-T,0) , 

provided suitable, actually optimal, regularity assumptions are made on the partial 
map {x,t) 1-^ a{t,x,z). Here C K" is an open subset, n > 2 and T > 0. 
Specifically, we shall consider a Caratheodory vector field a : fiy x R" — >■ M" which 
is C^-regular in the third variable and satisfying the following parabolicity and 
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continuity conditions: 

r |a(a;,t,z)| + |9.a(x,t,z)|(|z|2 + s2)i/2 <i(|z|2 + 52)(p-i)/2 

(1.4) I ,.{\z\^ + .s^)(P-^^/^\e<{dMxAz)^,0 

[ |a(a;,t,z)-a(xo,<,z)| <Lw(|x-a;o|)(|z|2 + s2)(P-i)/2 

whenever z,^ G M", x,a;o G fi, t G (— T, 0), where < < L arc positive num- 
bers and p > 2. The symbol uj(-) denotes a modulus of continuity meaning that 
uj: [0,oo) — > [0,1] is a nondecreasing concave function such that uj{0) = 0. In the 
following s > is a parameter that will be used to distinguish the degenerate case 
(s = 0), which covers the model equation in (|l.ip . from the non-degenerate one 
(s > 0); the analysis made in the following will see no difference between these two 
cases. In the rest of the paper we shall assume that the partial map 

a(x, t, z) 



(|z|2 + s2)(P-l)/2 

is Dini-continuos in the sense that 



(1.5) 



Jo Q 



This assumption is optimal for the estimates we are going to derive in the following; 
see the comments at the beginning of Section [1.31 We anyway remark that, every- 
where in this paper, we only assume measurability of the partial map t — a(x, i, z), 
in other words we assume that time coefficients are merely measurable. Yet, in this 
paper we shall always consider the case 

p>2 

as the case p < 2 has already been treated elsewhere [55] and involves an analysis 
which is different and somewhat simpler than the one which is necessary here. For 
more notation we refer the reader to Section 2 and to the rest of this introductory 
section. 

Remark 1.1 (On the notion of solution). Throughout the paper, when considering 
weak solutions to p.3p and unless otherwise stated, we shall mean energy weak 
solutions. An energy weak solution u belongs to the parabolic energy space, i.e. 

(1.6) u G C"(-r, 0; L'^in)) n LP{-T, O; W^'P{n)) , 

and it is a distributional solution to (|1.3p in the sense that 



(1.7) — / uipt dx dt + / (a{x,t, Du), Dip) dx dt =: / ipdfi 

holds whenever ip G C'^{^t)- In view of the available approximation theory we 
shall assume that /i G without loss of generality, while upon letting /i[Rn+i\ny = 
0, we shall finally consider the case 

(1.8) /iGLi(M"+i). 

Assumptions (|1.6p and p.Sp will be then finally removed in Section [TT^ There we 
shall deal with solutions to general measure data problems, where both (|1.6p and 
(|1.8p are not any longer in force. In other words, we purse the usual path that 
consists of first deriving a priori estimates for more regular problems and solutions, 
and then recovering the general case by approximation. Notice that under the 
assumptions p.6p and p.8p by standard density arguments the identity in (|1.7p 
remains valid whenever ip G W'^''p{Q.t) n L°°{Q.t) has a compact support. 
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1.1. Intrinsic and explicit Riesz potential estimates. Very recently, in [571 
[551 12H] for the case p > 2 and in [T3| for the subquadratic one, it has been shown 
that, surprisingly enough, tlie regularity theory of possibly degenerate quasilinear 
elliptic equations of the type 

— div a{Du) = /i 
completely reduces to that of standard Poisson equation 

(1.9) -Au^/i 

up to the C^-level, i.e. up to the gradient continuity. Moreover, in some sense the 
regularity theory can be actually linearized via Riesz potentials. In particular, the 
gradient of solutions can be pointwise bounded via classical Riesz potentials exactly 
as it happens for solutions to (jl.9|) . i.e., the inequality 

(1.10) \Du{xQ)y-^ <clY\xo,r) + c{i {\Du\ + s)d'j^ 

\JB{xo,r) J 

holds for a.c. point xq, where 

ir[xo,r) := / — 

Jo Q e 

denotes the standard truncated Riesz potential of Our aim is to build a related 
theory for general degenerate parabolic problems of the type in (jl.ip and (jl.3|) . The 
main challenge here is to match the anticipated a priori Riesz potential estimate 
with the inhomogeneous nature of equations such as (|l.ip : it will be indeed part 
of the work to find the proper formulation, suited to the geometry of the equations 
considered, making this possible. We also remark that, even when applied to the 
stationary case, our results turn out to be more general than those contained in 
[27l [28] since the equations considered here are also allowed to have coefficients. 
The ultimate outcome is that, once again, the (spatial) gradient regularity theory 
of solutions to (jl.3p can be unified in a natural way with the one of the usual heat 
equation (jl.2[) . The analysis here unavoidably involves the concept of the intrinsic 
geometry, introduced and widely employed by DiBencdetto [HI 1461 [2] . According to 
this principle, the lack of scaling (for p 7^ 2) of equations as 

(1.11) ut ~ div {\Du\P~^ Du) = 

can be locally rebalanced by performing the regularity analysis of the solution 
on certain special cylinders adapted to the solution itself, indeed called intrinsic 
parabolic cylinders. More precisely, instead of using standard parabolic cylinders 

(1.12) Qr{xo, to) B{xo, r) x [to - r^, to) , 

one uses cylinders whose time-length is stretched accordingly to the size of the gra- 
dient on the cylinder itself. In other words, one is lead to consider cylinders of the 
type 

(1.13) Q^{xo,to)-B{xo,r)xito-X''-Pr^,to), A>0, 
on which it simultaneously happens that a condition of the type 



(1.14) / 

JQ: 



\Du\ dxdt < A 

'Qi:{xo,tQ) 

is satisfied. The use of the word intrinsic stems from the very basic fact that 
the parameter A appears on both sides of (jl.l4p . Ultimately, this has the effect 
of rebalancing the local anisotropic character of the equation allowing for proving 
homogeneous regularity estimates: in some sense, the equation (jl.lip looks like 
the heat equation when considered on Q^{xo,to). For instance, when considering 
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standard parabolic cylinders, for solutions to it is only possible to prove 

bounds of the type 

(1.15) sup \Du\<c{n,p)j (\Du\ + s + ly^^ dx dt , 

whose lack of homogeneity precisely reflects that of the equation. In this sense 
the previous estimate is natural. When instead considering intrinsic cylinders with 
(|1.13p - p.l4p being in force, estimates become dimensionally homogeneous: 

/ \ i/(p-i) 

(1.16) c{n,p)if {\Du\ + s)P~^ dxdt] < \Du{xo,ta)\ < X. 

Both (|1.15p and (|1.16p are basic results of DiBenedetto [S] while we just remark 
that intrinsic geometries are nowadays a basic and common tool to treat degenerate 
parabolic equations pQ [THl HD] . 

The previous considerations, together with ()1.15[) - (jl.l6[) . are actually the starting 
point for proving the desired potential estimates. Let us see how. Beside the usual 
caloric (truncated) Riesz potentials built upon standard parabolic cylinders as in 
([TTT^ . that is 

(1.17) K{xo,to;r) := , < P < N := n + 2 , 

Jo Q Q 

we introduce the intrinsic Riesz potentials as 

r |Ai|(Q^(a;o,to)) dg 



(1.18) ^{x^M;r) := / 

Jo 



g^-P g 



where N is called, as usual, the parabolic dimension. Note that in such a way we 
have l^{xQ,tQ]r) = 1^ ^{xQ,tQ]r). At this stage the word "intrinsic" merely refers 
to the fact that the additional parameter A has been considered in the definition 
in (|1.18p . while at the moment no local linkage with solutions of the type in (jl.l4p 
has been considered yet. This will come in a few moments: indeed, the right way 
to give an intrinsic formulation of the linear potential bounds is inspired by (jl.l6p 
and it is given in the following: 

Theorem 1.1 (Intrinsic Riesz potential bound). Let u he a solution to (jl.3p under 
the assumptions (jl.4p - (jl.5p and (|1.8p . There exist a constant c > 1 and a radius 
Rq > 0, both depending only on n,p,iy,L,Lu{-), such that the following implication 
holds: 



(1.19) 




cl^^(a;o,io;0 +c f {\Du\ + sf''^ dx dt 




\Du{xq, to) I < A 



whenever Q^(^Oiio) C Ht, {xa,to) is Lebesgue point of Du, and r < Rq. When the 
vector field a(-) is independent of x, no restriction occurs on r, i.e., R^) — oo. 

Note that, as in a sense it was a priori required, (|1.19p allows to recover (|1.16p 
when /i = 0; this is a first sign of the fact that (jl.l9p is the "correct intrinsic 
extension" of (|1.10p . As a matter of fact Theorem 11.11 implies a gradient linear 
potential estimate involving standard Riesz potentials. Surprisingly enough, this is 
of the same type as the one which holds for the standard heat equation; moreover, 
when jU = 0, this reduces (|1.15l) . We indeed have the following: 



RIESZ POTENTIALS AND NONLINEAR PARABOLIC EQUATIONS 



5 



Theorem 1.2 (Riesz potential bound in classic form). Let u be a solution to (|1.3p 
in Qt under the assumptions (|1.4p - (jl.5[) and (|1.8p . There exists a constant c, 
depending only on n,p, z^, X, w(-), such that 

\Du{xo,to)\ <cl[^^{xo,to;r)+c-f {\Du\ + s + If'^ dx dt 

J Qr{xQ,to) 

holds whenever (xQ,to) G ^^t is o, Lebesgue point of Du and whenever QrixQ,to) C 
fix is a standard parabolic cylinder such that r < Rq; here Rq is the radius intro- 
duced in Theorem ] 1. 11 

For the case p = 2 the previous estimate recovers the main parabolic result in 
[l2] . An immediate consequence of Theorem 11.21 is the following global bound via 
classical, non-truncated caloric Riesz potentials: 

Corollary 1.1. Let u he a weak solution to the equation 

(1.20) ut - diva{t,Du) = fj. in ]R"+^ 

under the assumptions (jl.4p ; moreover, assume that the global integrability u S 
LP~^{—oo,tf)]W^'P~^{W^)) holds for to G M. There exists a constant c, depending 
only on n,p,iy,L, such that the upper bound 

J{t<to} dpi,r{[x,t), [Xo^to))^" 

holds whenever (xo,io) ci Lebesgue point of Du. 

Remark 1.2. In Theorem 11.11 dp=r(-) denotes the standard parabolic distance in 
M"+^, which is defined by 

rfpar((a;,i), {xQ,to)) := max||a; - xq\, ^/\t - <o|| • 

The previous result shows that Theorems 11.11 and 11.21 play the role of the usual 
representation formulae via heat kernels for solutions to the heat equation. In recent 
years there has been a large activity devoted to the understanding the extent to 
which heat kernel estimates are still valid when passing to more general settings, 
as for instance Lie groups and manifolds [TTJ ST] . In this paper we are interested, 
in a dual but yet related way, to see the extent to which estimates as those implied 
by wcU-bchaving heat kernels can be recovered in the nonlinear degenerate setting. 
Our results also connects to a recent fact observed in |34j . and concerning the 
p-superharmonicity of linear Riesz potentials. 

The proof of Theorem 11.11 opens the way to an optimal continuity criterion for 
the gradient still involving only classical Riesz potentials and that, as such, is again 
independent of p. 

Theorem 1.3 (Gradient continuity via linear potentials). Let u be a solution to 
(II. 3p in flT under the assumptions (|1.4l) - (jl.5p and (|1.8p . // 

lim I't{x,t:r) = 

locally uniformly w.r.t. {x,t), then Du is continuous in fix- 

An important corollary involves Lorentz spaces: 

Corollary 1.2 (Lorentz spaces criterion). Let u be a solution to (|1.3p in Qt under 
the assumptions p.4p - (|1.5p . Lf ji ^ L{N, 1), that is if 

/•oo 

/ \{{x,t) e nr : Hx,t)\ > A}|i/^dA < oo , 
Jo 

then Du is continuous in fix- 
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Corollary 11.21 substantially improves the ones in [53] claiming only the bound- 
edness of the gradient under the assumption /i e L{N, 1); see also [B] for a related 
global boundedness result in the elliptic case. It might be interesting to note how 
the above result naturally extends to the parabolic case the classical gradient conti- 
nuity results valid in the elliptic case, starting from those available for the Poisson 
equation — Au = /i in domain of M". For this it is known that the condition 
11 S L(n, 1) is a sufficient one for the gradient continuity. This is in turn related to, 
and indeed implied by, a classical result of Stein [55] that claims the continuity of a 
function / whenever its distributional derivatives belong to L(n, 1). CoroUarv 11.21 
gives the precise nonlinear parabolic analog of such classical facts. As expected, the 
space dimension n is replaced by the parabolic one iV = n + 2, which is naturally 
associated to the standard parabolic metric. 

Preliminary to the proof of the continuity criterion, there is another result which 
claims the VMO gradient regularity under weaker assumptions on the measure /i. 

Theorem 1.4 (VMO gradient regularity). Let u he a solution to p.3p under the 
assumptions (|1.4p - (|1.5p and p.Sp . Ifli{x,t;r) is locally bounded in Qt md if 

(1.21) iin,M%(^ = 

locally uniformly in w.r.t. {x,t), then Du is locally VMO in fix, that is 



lini sup -f \Du — {Du)q^ \ dxdt ~ 

R^Q r<R.Q,.<ZQ' J 



for every open subset Q' (e fix- 

1.2. General measure data problems. Solutions to measure data problems are 
usually found by approximation procedures via solutions to more regular problems. 
These are of the type 

(1.22) {uh)t - div a{x,t, Duh) = fih^ , heN, 

where (jl.6|) holds for Uh and fih is usually a convolution of fi. The point is that 
solutions to measure data problems do not belong, in general, to the energy space. 
This section is also aimed at justifying that we may actually work under the appar- 
ently additional assumption p.6p . More precisely, the exact definition of SOLA, is 
given in the following: 

Definition 1 ( [H dU [52] ) . A SOLA (Solution Obtained as a Limit of Approxima- 
tions) to ([rg is a distributional solution u e LP-^{-T,0;W^-P-^{n)) to ([O]) in 
VIt, such that u is the limit of solutions 

Uh e C"{~T, 0; L'^in)) n LP{~T, 0; W^'P{n)) 

to equations as (fr22]) in the sense that uh u in LP~^{-T,Q] W^^p-'^(SI)), L°° 3 
fih fJ- weakly in the sense of measures and such that 

(1.23) limSUp \fih\{Q) < iMl(LQJpar) 

h 

for every cylinder Q = B x {ti,t2) ^ f^T, where i? C is a boimded open subset. 

In the right hand side of (|1.23p appears the symbol [QJpar, which denotes the 
parabolic closure of Q defined in (|2.ip below. For more on this kind of solutions see 
Rcmark l5.1l b clow : in particular, requiring (jl.23p is neither unnatural nor restrictive. 
Our estimates remain valid for SOLA and, in fact, the following holds: 



Theorem 1.5. The statements of Theorems M.li uTB. \1.3\ and \1.4\ continue to hold 
for SOLA u e LP-^{-T,0;W^^P-^{n)) to ([TS]), under the only assumptions ([13)- 
(|1.5p . As a consequence, the results in Corollaries \1.1\ and \L3\ hold for SOLA as 
well. 
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1.3. Comparison with nonlinear estimates. Theorem 11.11 improves the previ- 
ously potential estimates via nonlinear potentials |25j , bringing them to the desired 
optimal level. Based on elementary dimension analysis we conjecture that the result 
of Theorem 11.11 cannot be improved by the use of any other nonlinear potential. 
Theorem 11.11 is optimal also with respect to the regularity assumed on the coeffi- 
cients dependence x a{x, •), that is (|1.5|) . Indeed, already in the linear elliptic 
case 

div {d,{x)Du) = , 

Dini-continuity of elliptic coefficients matrix a{x) is essential in order to get gradient 
boundedness. Merely continuous coefficients are not sufficient to ensure that the 
gradient belongs even to BMO, see [TB] . 

Now, let us see how Theorem 11.11 improves the previously known estimates via 
nonlinear potentials. In [35] a Wolff potential type gradient bound has been ob- 
tained for equations without coefficients of the type Ut — div a{Du) = fi. More 
precisely, in [25] we introduced the following intrinsic Wolff potentials: 



^l\iQ^g{xo,to)) 

X2-pgN-l 



l/(p-l) 



dg 

— , N ^ n 
Q 



A > 0. 



See [21 [T5J dH] for more on Wolff potentials. We then proved the existence of a 
universal constant Cw = Cu,(7i,p, L) for which 



(1.24) 



c„,W^(a;o,io;0 



i/(p-i) 



Q^(3:o,to) 



{\Du\ + sy-^ dxdt 



< A 



=^ \Du{xo,to)\ < A 
holds. Let us now show that (|1.19l) implies (|1.24p . for 

c.:= [4^+i(log2)2-Pc]'/^^-'^ 



(1.25) 



4^c 



and radii r < Rq, where c and i?o are the constants appearing in the statement of 
Theorem ll.il (no restriction on radii appears in the case of equations as in (jl.20p V 
With ri ~ r/2* for integers i > 0, Holder's inequality for series (as it is p > 2 here), 
gives 



E 

i=i •^'^•+1 



< 



< 



/i|(g^(a:o,to)) dg 

oo i..i/^A 



lMl((3n(a;o,io)) 



log 2 



E 



\fi\{Qi:^{xo,to)) 

„N-1 



1/(P-1)' 



1=1 



< 4"(log2)2-P 



E 



^^\iQ^{xo,to)) 

nN-l 



P-1 



1/(P-1) 



dg 

g 



p-1 



4^(Iog2)2-fA2-''[W^(zo,to;0] 



p-1 



< 4^(log2)2-Pci,~?'A < 



2c 



Notice that to derive the second-last estimate we have use the inequality in the first 
line of (|1.24p . while in the last estimate we have used (|1.25p . Using the inequality 
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in the last display, again (|1.25p and finally the left hand hand side of (|1.24p . a 
standard manipulation gives 

cI';^{xo,to;r/2) + ci f {\Du\ + s)p-Ux dt] <X 

SO that the right hand side inequality in (|1.24p . that is \Du{xo,to)\ < X, follows 
applying Theorem ll.il 

The improvement from (|1.24p to (|1.19p is rather strong both from the viewpoint 
of the theoretical significance - as now the regularity theory of quasilinear equations 
is unified with that of the heat equation up to spatial gradient continuity - and 
from the one of the consequences. For instance, when looking for sharp criteria for 
establishing Du € L°° (and eventually in C°) in terms of Lorentz spaces, the result 
in (fLM)) gives that 

(1.26) ii<EL{N,ll{p-l))=^ DueL°° , 
where we recall that 

/i e L(iV, l/(p- 1)) / Al^|{(x,t) e : |/i(a;,i)| > A}|W^rfA < c». 

Jo 

The condition of Corollarv ll.2l is clearly stronger that the one in (ll.26p . as L{N^ 1) C 
L{N,l/{p — 1)), this inclusion being strict for p > 2. For further properties of 
Lorentz spaces we refer for instance to [H]. Moreover, it is easy too see that more 
refined criteria in terms of density /concentration are provided by (|1.19p with respect 
to (|1.24p when fj, is genuinely a measure. We also remark that Wolff potentials play 
a major role in the analysis of the fine properties of quasilinear equations (see 
for instance [T71 [T8j [39l 001 145]); since the estimates contained in this paper are 
stronger than those involving Wolff potentials, we expect they will have a similar, 
if not stronger, impact in the future. 

1.4. Techniques. Finally, a few comments on the methods used in this paper. 
Several new ingredients are needed to deal with the parabolic case with respect 
to the previous elliptic one l24|, and the proofs depart considerably from those 
proposed before. The proof of Theorem 11.11 involves a very delicate, double-step 
induction procedure based on a few ingredients that re-shuffle, in a pointwisc man- 
ner, some classical methods used in linear Caldcron-Zygmund theory and combine 
them with the use of intrinsic geometry. Extensive use of nonlinear potential the- 
oretic methods and regularity theory is made throughout. Let us briefly describe 
the heuristic used here by specializing for simplicity to the model case (jl.ip and 
considering yu G L^; the essence relies in careful procedure that allows to "linearize 
the equation" and control the possible degeneracy in a precisely quantified way at 
every scale. The following argument will be purely formal. We consider a dyadic 
shirking sequence of intrinsic cylinders for i > 

•■•Qr.+i(a;o,io) C Q^^(.To,to) C Q^^_^{xo,to) . . . n := a\ 

where u e (0, 1) is a constant depending only on n^p and A is as in (jl.l9p . A suitable 
exit time argument, together with very careful regularity estimates for solutions to 
homogeneous equations, gives 

(1.27) A — "quantified error" < I Dw I on Q^. (xq, to) for i large enough . 

This is something that in a way we can always assume, starting from an exit time 
index, otherwise we are going to get an opposite inequality for the integral averages 
(^o,to)l ^ that eventually leads to an immediate proof of \Du(xQ,tQ)\ < 
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A, and therefore of (|1.19p . Assuming (|1.27p leads to implement a delicate iteration 
procedure whose finally outcome is the following inequality: 

(1.28) -f \Du\P-Ux dt < XP^^ 

that again implies (|1.19p . Note that proving (|1.28p not only allows to implement 
the iteration but also allows to use, at each scale, the intrinsic geometry methods 
(compare with (jl.l4p '). As emphasized, the key to the proof of Theorem 1 1.1 1 is the 
lower bound in (jl.27p ; let us now give a formal but yet convincing argument showing 
how a condition as (jl.27p allows to get (|1.19p and why intrinsic Riesz potentials 
and and conditions as (|1.19p naturally occur. Let us consider then (jl.27p to be 
satisfied on (5^(xo,to) with a null error, i.e., A < \Du\^ and let us assume the first 
inequality in (|1.19p . The lower bound A < \Du\ in turn allows to gain coercivity 
enough to treat the equation in display (jl.ip as a heat equation with a coefficient, 
that is ut — div {Xp~^Du) ~ fi that we can rewrite as 

X^-Put ~Au = X^-Pfx in Q^{xo, to) . 

Now the effect of the use of intrinsic geometry and of the intrinsic Riesz potential 
shows-up. Changing variables and introducing 

/ ,^ u{xo + rx,to + X'^-PrH) ^ 2-p 2.^ 

v[x,t) -.^ , fi[x,t) := X Pri^iixo + rx,tQ + X ' rt) 

r 

for {x,t) e Qi = Bi X (—1,0), we have 

(1.29) Vt-Av = fl. 

Next, we apply the standard Riesz potential bound for solutions to (|1.29p . that is 

/ X i/(p-i) 

(1.30) |i:»w(0,0)|<I^'^'(0,0;l)+ -/ \Dv\P-Uxdt] 

\JQi(o.o) J 

Changing variables back to fi wc notice 

lf'(0,0;l) = X^-Pr [ / \fi{xo+rx,ta + X^-Pr^t)\dxdtdQ 

■JO JQeiOfi) 

= X^^Pr f -f \fi{x,t)\dxdtdg 

= X^'P [ -f \n{x,t)\dxdtdg^X^-PI'^^{xo,to;r)<X, 

Jo JQ^ixo,to) 

where in the last inequality we have used the first line in (|1.19p . Finally, scaling 
back to u, using the previous inequality together with the first line of (|1.19p . we 
conclude with \Du{xo,to)\ = |-D?;(0,0)| < A, that is the proof of (|1.19p . The one 
outlined in the last lines is only a heuristic argument used to show how intrinsic 
Riesz potentials play a decisive and natural role in this context, but its rigorous 
implementation is highly nontrivial and involves a refined double induction argu- 
ment that exploits rather subtle aspects of regularity theory of degenerate parabolic 
equations. Several tools are used here. One of the main points is that the anal- 
ysis of the relevant iterating quantities must be performed at two different levels, 
using different energy spaces. Indeed, since we are dealing essentially with mea- 
sure data problems, the natural spaces involved are larger than LP{—T,0;W^'P). 
This, together with the lack of reverse Holder type inequalities and homogeneous 
estimates which is typical when dealing with degenerate parabolic equations, re- 
flects in a simultaneous use of two different spaces, namely L^{—T,0;W^-^) and 
LP~^{—T, 0; W^'P^^). Eventually, a very delicate interplay between local regularity 
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of solutions to homogeneous equations and comparison estimates must be exploited 
in the framework of intrinsic geometries thanks to exit time arguments and the use 
of intrinsic Riesz potentials. The proof of Theorem 11.11 eventually opens the way to 
the continuity analysis and in particular to Theorem II. 31 For this we shall readapt 
the iteration procedure of Theorem II. II to estimate oscillations rather than the size 
of the gradient. This in turn imposes to consider a priori infinite many exit times 
arguments used to control the degeneracy of the equation via the oscillations of the 
gradient, and vice- versa. 

2. Preparations 

2.1. General notation. In what follows we denote by c a general positive con- 
stant, possibly varying from line to line; special occurrences will be denoted by 
ci , C2 , ci , C2 or the like. All these constants will always be larger or equal than one; 
moreover relevant dependencies on parameters will be emphasized using parenthe- 
ses, i.e., ci = ci{n,p,v^ L) means that ci depends only on n,p^iy,L. We denote by 
B(xo,r) := {x G R" : |a; — a;o| < ?'} the open ball with center xq and radius r > 0; 
when not important, or clear from the context, we shall omit denoting the center as 
follows; Br = B{xo, r). Unless otherwise stated, different balls in the same context 
will have the same center. We shall also denote B = Bi = B{0, 1) if not differently 
specified. In a similar fashion standard and intrinsic parabolic cylinders with vertex 
(xq, <o) and width r > have been defined in p.l2p and (|1.13p . respectively. When 
the vertex will not be important in the context or it will be clear that all the cylin- 
ders occurring in a proof will share the same vertex, we shall omit to indicate it, 
simply denoting Qr and for the cylinders in (|1.12p and p.l3p . respectively. We 
recall that if Q = Ax {ti,t2) is a cylindrical domain, the usual parabolic boundary 
of Q is 9parQ := {Ax {ti})\j{dAx [^1,^2)), and this is nothing else but the standard 
topological boundary without the upper cap ^ x {^2}- Accordingly, we define the 
parabolic closure of Q as 

(2.1) LQJpar :=QUaparg. 

With O C R""*"^ being a measurable subset with positive measure, and with g: O ^ 
M" being a measurable map, we shall denote by 



(5) 



-f g{x,t) dx dt :~ -—- j g{x,t)dxdt 
Jo 1^1 Jo 



its integral average; here \0\ denotes the Lebesgue measure of C. A similar notation 
is adopted if the integral is only in space or time. In the rest of the paper we shall 
use several times the following elementary property of integral averages: 

(2.2) \g-{9)o\''dxdtj <2\J- \g--i\'idxdt 

whenever 7 G M" and q>l. The oscillation of g on O is instead defined as 
osc g := sup \g{x,t) - g{x,i)\ . 

^ {x,t),{x,i)€0 

Finally, we remark that we shall denote the partial derivative with respect to time 
of a function u both by ut and by dtu; moreover, the letter A will always denote a 
positive number. Further relevant notation is at the beginning of the next section. 

2.2. Comparison maps. The basic setup in this section is tailored to the needs 
of the proof of Theorem 11.11 and subsequent results. Therefore we shall consider u 
to be an energy solution to (|1.3p under the assumptions (|1.4p - (|1.5p and (|1.8|) until 
the end Section only in Section [^751 we shall discuss the general case, thereby 
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treating SOLA and discarding assumption (|1.8p . With a point {xo,to) G i^r being 
fixed, and given an intrinsic cylinder of the type 

QUxo,to) = B{xo,r) x [to - X''-Pr^,to) 

such that Q2r{xo,to) C fix, we consider a family of nested parabolic cylinders 

(2.3) Qj = Bj X T, = B{xo, r,) x {t^ - \^~^r]M) ^ , := , 

for a fixed decay parameter cr G (0, 1/4); notice that we always have 

■■■Qj C ^Qj-i C Qj_i . . . 
Accordingly, we consider their dyadic, parabolic dilations 

tQj = Q^,^(xo,to) = tB, X tT, = B{xQ,Trj) x (to - A^-f (rrj)', to) 
for T > 0; notice that here, slightly abusing the notation, we are denoting 

A similar notation will occur several times in rest of the paper. Now, let 

w, e C^{T,-L\B,)) n LP{T,-W^^P{B,)) 
be the unique solution to the Cauchy-Dirichlet problem 

^2 I dtWj - div a(x, t, Diu^ ) = in Q^- 

I Wj ^ u on dpaiQj . 

After having defined Wj , we also define 

v,eC"{^T,;L'{^B,))nLP{^T,;W'-P{^B,)) 

as the unique solution to the frozen Cauchy-Dirichlet problem 

{dfVj — div a{xo,t, Dvj) =0 in ^Qj 
Vj = Wj on apar {\Q])- 



(2.5) 



2.3. A priori estimates for comparison maps. We now derive various a priori 
estimates for Wj and Wj, starting from L°°-bounds. When turning our attention to 
Wj we need to use the results recently established in [5^, that allow to deal with 
equations with non-constant, Dini-continuous spatial coefficients. We start with a 
statement in terms of intrinsic geometry. 

Theorem 2.1 (Intrinsic gradient bound). Let Wj he as in (|2.4p . There exists a 
positive radius Ri = Ri{n,p,v, L,uj{^)) and a constant ci = ci(n,p, i^, L) such that 
if g ^ (0, i?i) and 

Q]°{xiM) ■■^B{xi,e) X (ii - Ao^Pp^ ii) c Qj , 

then the implication 

\ i/p 

\Dwj\+ sf dxdt] < Ao \Dwj{xi,ti)\<\Q 
holds. 

Proof. This result has been proved and used in [291 Theorem 1.1, Theorem 1.3, 
Theorem 4.1]. The proof is exactly the one given in the proof of [29j Theorem 1.1], 
once [121 Lemma 4.3] is used instead of [IHl Lemma 4.2], see also [IHl Remark 4.1]. 
We also remark that no restriction on g is necessary when the vector field a(-) is 
independent of x. □ 
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The pointwise bound of Theorem l2.1l can be turned mto an L°°-bound of exactly 
the same type proved by DiBenedetto |2] for equations with no coefficients, see also 
Theorem l2. 21 below. Since we are going to cover the case of equations with measure 
data, where solutions with low degree of integrability naturally appear, we need 
to lower the p-integrability exponent to (p — 1) to get the correct form of a priori 
estimates. All this is done in the next 

Corollary 2.1. Let Wj he as in (j2.4p and Ri as in Theorem \2.1[ There exists a 
constant C2 = C2(n,p,i', L) such that if r G (0,i?i), then 

sup \Dw, I + s < C2(A + ,5) + ^^^^,^f^ .„^, / {\Dwj I + sy-^ dx dt 

holds whenever t„i G (0, 1). In particular, we have 

sup \Dwj\ + s < C2(A + .s) +C2A^"P f {\Dwj\ + sf'^ dxdt. 

Proof. Define 

1 2r^ r 
Ao = Ao(t, t') - sup \Dw, \+\ + s + 'jr;^^^-^ f i\Dw,\ + sY"' dx dt 

whenever r„j < t' < r < 1, where ci = ci{n,p,i>, L) is as in Theorem 12.11 As 
Ao > A and p >2, we clearly have that for 5 :~ t — t' the inclusion 

Qt,{xiM) ■■= B{xuSr,) X (ii - Xl-P{Srjf,ti) C rQ, = tQ,^^ (x-q, to) 

holds whenever G r'Qj. Furthermore, by using the very definition of Ao we 

may estimate 

ci I f {\Dwj\+ s)P dxdt] 

\JQi°.{x^,t^) J 



/ ilDwjl+sf^dxdt] 

<c,2'/PX'JP ( ^ ( I {\DwA+sr-Uxdt] 

ij {\Dwj\ + s)P-^dxdtj 




Therefore Theorem 12.11 implies that \Dwj{xi,ti)\ < Aq. But this holds for all 
{xi,ti) G t'Qj and thus 

sup \Dwj I < - sup \Dwj \ + X + s+ ^ -[■ {\Dwj \ + s)''"^ dx dt 

T'Qj ^ TQj [T - T )"+ Jq. 

follows. Lemma 12.11 below applied with lp{t) = sup^g^. \Dwj\ then concludes the 
proof by properly choosing the constant C2 . □ 

The nest result is a classical iteration lemma for the of which we refer to [?n 
Lemma 6.1]. 
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Lemma 2.1. Let tp : [t„i, 1] — > [0, oo), with t„j G (0, 1), be a function such that 
</3(r') < -<^(t) + K + ——r holds for every r„ < t' < t < 1 , 

where B, K > 0. Then (^(r„) < c{n)K + (1 - r„0"^"+^^S. 

Corollary 12.11 obviously holds for Vj too, and in this case it is a by now classical 
estimate of DiBenedetto [9] , as already mentioned above. See also for example [H 
[25] for similar bounds. We report the statement for completeness. 

Theorem 2.2. Let Vj be as in (|2.5|) . For a constant C3 = cy,[n^p, v^L) we have 
sup \Dv.j\ + s< C3(A + s) + caA^^P / {\Dvj\ + sf^^ dxdt. 

We now pass to give oscillation estimates for Wj and Vj . The next result provides 
a gradient oscillations estimate for solutions to homogeneous equations with Dini- 
continuous coefhcients. 

Theorem 2.3 (Continuity estimate). Let Wj be as in (|2.4p . then Dwj is continu- 
ous. Moreover, assume that 

(2.6) sup \Dwj\+ s < AX 

holds for some A > 1. Then, for any S G (0,1) there exists a positive constant 
ui = (Ti{n,p, V, L, A, S, uj{-)) G (0, 1/4) such that 

(2.7) osc Dwj < SX . 

Proof. The starting point of the proof is the work in [52] , where the continuity of the 
gradient of solutions to parabolic equations as in (|1.3p has been proved under the 
assumption of Dini-continuity of the space coefficients; this by the way immediately 
implies the continuity of Dwj claimed in the statement. What we need here is a 
quantitative bound on the oscillations of Dwj . To this aim, let us briefly recall the 
main arguments in [551 proof of Theorem 1.3], where the continuity properties of 
Du are formulated and proved in terms of the auxiliary vector field 

V{z) :=(|z|2 + s2)(P-2)/4^ 

and the related field V{pwj). For the use of such maps in the present context we 
refer to [52]; the only property we shall use here is the following inequality: 

W -1^.-. \V{zi)-V{z2)\ 

that holds for c„ = Cy{n,p) and for all vectors ^1,2:2 G K" which are not simul- 
taneously null; see for instance [37) . By following the arguments developed for 
[5^ (5.15)] it can be proved that for every e G (0, 1) there exists a positive radius 
Re = Re{n,p,v,L,ijj{-),£) G (0,1/16) such that 

(2.9) \(y{Dw,))Q.>.^^:,^ - {V{Dw,))q.>.^^;,^\ < {AXf/'e 
and 

(2.10) (-f \V{Dw,)-{V{Dw,))Q..^.,^i^\^dxdt] <iAXr/^e 

hold whenever (i,<) G jQj and < t < g < Rerj] notice that R^ is in particular 
independent of A, A and the considered cylinder Qj. Letting t — > in (|2.9p and 
recalling that V{Dwj) is continuous yields 

(2.11) \V{Dw,{i,i)) - iV{Dw,))QA.^^^i)\ < {AXy/'e V g G {0,R,r,] . 
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We are now ready to finish the proof with the ehoice 



c„2P/2-i48^Ap/2' 32 ■ 

The constant Cy is the one appearing in (|2.8p Notice that the dependence of tri 
upon n,p, I/, L, j4, (5, aj(-), as described in the statement, appears through the one 
imphcitly contained in i?^. Take now (ij,s),{x,t) G uiQj\ we can assume that t>s 
otherwise we can exchange the role between the two points in the next lines. It 
obviously follows that 

Using this last fact, thanks to Jensen's inequality, the one in display (|2.10p and 
using also (|2.2p . we have 



/ \V{Dw,)~[V[Dw,))qa. ,^,i)\dxdt 

2-f \ViDw,)-{ViDw,))QA. (,^i)\dxdt 

/ \V{Dw,)-{V{Dw,))qa>. ^,~t)\dxdt 



< 16^ 



< 16^ 

(i \V{Dw,)-{V{Dw,))qa. 

< 16^(AA)f/2e. 

By using the previous estimate and ()2.1ip (actually used also for {y, s) instead of 
(x,i)) together with triangle inequality we easily gain 

(2.12) \V{Dw^{S;,i)) ^V{Dwjiy,S))\ < iS'^ iA\)P/^e . 
We are now ready to show (|2.7p proving that 

(2.13) \Dwj{x,i) - Dwj{y,S)\ < SX 

whenever (y, s),(i,i) G oiQj. First of all, observe that we can assume that ei- 
ther \Dwj{x,t)\ > SX/2 or \Dwj{y, s)\ > dX/2 holds otherwise we are done. The 
inequalities in (|2.8p and (|2.12p then imply (|2.13p as follows: 

\Dwj{x,i) - Dwjiy,S)\ < c„48^(^A)P/2(JA/2)i-P/2^ = SX 

and the proof is complete. □ 

Remark 2.1. The proof of the previous result allows in fact, together with the 
argument given in |29| , to get an explicit modulus of continuity for the gradient of 
solutions of equations with Dini-continuous coefficients. Indeed, the choice of the 
radius Refj making (|2.9p is made to meet a condition of the form 

r^^'^^ do 
u:{R,rj)+ / w(^^) — < - 
Jo g c 

for some positive constants b and c depending only on n^p^v, L, A. This gives a 
modulus of continuity involving a power of the function 

dg 



r t-^ aj(r) + / i^{g) 

Jo Q 

which is in accordance to the known results in the classical elliptic regularity theory. 
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We now collect a few results from [551 US] in order to provide oscillation estimates 
for the functions Vj. The next statement is a slight variant of Theorem 3.1]. 

Theorem 2.4. Let Vj be as in (j2.5p . Consider numbers 

A,B>1 and e G (0,1) . 

Then there exists a constant (72 G (0, 1/4) depending only on n,p,iy, L, A, B, e such 
that if 

(2.14) — < sup |Dwj | < s + sup \Dv^ \ < AX 
holds, then 

(2.15) / \Dvj - (-Dfj)rQ, \dxdt<ef \Dvj - {Dvj)i \ dxdt 
holds too, whenever r G (0, (72] ■ 

Remark 2.2. The essence of the previous result lies in the fact that once the 
bounds (|2.14p are satisfied, then solutions to evolutionary p-Laplacean type equa- 
tions satisfy elliptic type decay estimates as in (|2.17p when framed in the proper 
intrinsic geometry dictated by (|2.14p . Indeed, let us denote by E{f,Q) the usual 
excess functional 



(2.16) E{f,Q):^ -f \f-{f)Q 

JQ 



dx dt 



which is defined whenever / is an integrable function and Q a measurable set with 
positive measure; this functional gives an integral measure of the oscillations of / 
in a subset Q. Estimate (j2.15p now reads as 

(2.17) E{Dvj , rQj ) < eE{Dvj , i ) . 

Theorem l2. 41 gives the natural analog, when passing to the framework of degenerate 
parabolic equations of p-Laplacean type, of the classical results known for solutions 
to the heat equations. Indeed, Theorem 12.41 holds without assuming (|2.14p for 
solutions to p. 21) . This is a classical result of Campanato [5]. 



Using Theorem 12.41 it is possible to give a proof of the Holder continuity of the 
gradient of solutions to frozen equations, as for instance shown in |26[ Theorem 
3.2]; see also [ISJ Theorem 3.2]. 

Theorem 2.5. Let Vj be as in (j2.5p . For every A > 1 there exist constants C4 = 
C4{n,p, V, L, A) and a = a(n,p, u, L, A) such that 

sup \Dvj \+s<AX =^ osc Dvj < C4r"A Vt G (0, 1/4) . 

Remark 2.3. Theorems 12.4)15751 have been presented in [551 [55] actually for solu- 
tions to equations of the type 

(2.18) ut - div a{Du) = . 

On the other hand, by following the arguments in [55] [5B] it is not difficult to see 
that all the proofs carry out for solutions to equations of the type in display (|1.20p . 
Therefore Theorems 12.4112.51 apply to the functions Vj as well, that indeed solve 
equations as in (|1.20[) . 
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2.4. Comparison estimates. We start this section by a reformulation of a result 
established in [551 Lemma 4.1] and [52] (4.5), (4.6)]. We remark that the result there 
was presented only for equations without coefficients as in ([2TT8)) . Nevertheless, the 
proof works directly for general equations with merely measurable coefficients; the 
crucial point is the strict monotonicity in the gradient variable. 

Lemma 2.2. Let u he as in Theorem \1.1\ and Wj as in (j2.4p with j > 0. Let 
£ G (0, l/(n + 1)]. There exist constants Ci = Ci(n,p, i/, e) and C2 = C2{n,p,i>) such 
that 



(2.19) 



1/9 



\Du- DwA'^dxdt 



< ciA 



1 WQi 



(n+2)/[(p-l)n+p] 



holds for any 0<(7<]3— l + l/(n+l) — £. Moreover, the inequalities 



(2.20) 
and 
(2.21) 



sup 



\u — Wj \ dx < \fJ-\{Qj) 



/ 

Jq 



Du\ + \Dw,\)P-^\Du - Dwi\^ 

' ' -"^ ' -^dxdt<C2 



hold for any a > and ^ > 1 . 



We here recall a parabolic Sobolev-Poincare inequality that will be useful in the 
sequel; we refer to [5] Chapter 1, Proposition 3.1] for the proof. 

Proposition 2.1. Let v e L°°{Tj; L"''{Bj)) n L'^^ {Tj-W^'''^Bj)) for q2,m > 1. 
There exists a constant c depending only on n,q2,m such that the following in- 
equality holds for qi ~ q2{n + m)/n: 

\ 92/" 



dxdt < c 



Qj 



/ \Dv\'i'dxdt \ [sup / |w(x,- 



■)rdx 



Using the previous result and Lemma 12.21 we get another comparison estimate. 



Lemma 2.3. Let u be as in Theorem and Wj^i,Wj as in (|2.4[1 . with j > 1. 
Then, for any e e (0, l/(ri + 1)], these exists a constant C3 = C3(n,p, i^, e, cr) such 
that the inequality 



(2.22) 



1/9 



\u — I' dx dt 



< C3r,_iA 



1 ImKQj-i) 



-1 {n+p)/[{p-l)n+p 



'i-1 



holds whenever 0<q<p— 1+ p/n — £ and 

\ 1/9 



(2.23) 



iDwj^i — DwjV^ dxdt 



< C3A 



1 MQ. 



'j-i 



(ri+2)/[(p-l)n+p] 



holds whenever 0<q<p— l + l/{n+l)— £. 

Proof. To prove (|2.22p . we use Proposition 12.11 with the choice to = 1, 

p ^ n _ 1 77,e 



91 



p-l 



92 



-qi=p-l 



n n + 1 n + lrt + 1 

and u = M - Wj e L°°{Tj;L^{Bj)) n L'i^{Tf,W^''^^Bj)); recaU ([5:50)) . This yields 

1/91 
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l/n\ n/{n+l) 



< C 



j 





l/«2 




dx dt 




sup / \u — Wj\dx 









Let ci be as in Lemma be the constant corresponding to the choice en/{n + 1) 
instead of e. Substituting (|2.19p and (j2.20p into the previous estimate leads to 

\ 1/91 

\u-Wj\''^ dxdt] 



;"/("+!) 



< cc 



< C3r,_iA 



1 ImKQj) 



-| {n+2)/[{p-l)n+p] 



n/{n+l) 



l/n 



1 WQj-i) 



'3-1 



("+P)/[(P-1)"+P] 



which, together with Holder's inequality, proves (j2.22p . Here we have also used that 
— 1 = n + 1 and the identity 



n + 1 



n + 2 



1 



(p — l)n + p n 



n+p 



{p — l)n + p 



As for (|2.23p we instead argue as follows: 



1/92 



\Dwj^i - Di 



dx dt 



< 



\Q. 



j-ii 



\Qj 



f \Du- Dwjl'^^ dxdt 



1/92 



Du - Dwj^il''^ dxdt 

1/92 



< Ci 



< C3A 



IQ.I 



1/92 ^w-i 

, 3-1 



„w-i 



1 ImKQj-i) 



-| (n+2)/[(p-l)n+p] 



'j-1 



1 ImI(Q,-i) 



'j-l 



(n+2)/[(p-l),i+p] 



where we have repeatedly applied (|2.19p and used the fact that we arc assuming 
p > 2. Now (|2.23p follows, again by Holder's inequality as e is arbitrary. □ 

The following lemma provides one of the key estimates to obtain Theorem 11.11 

Lemma 2.4. Let u be as in Theorem \1.1\ and Wj^i,Wj as in (|2.4p . with j > I. 
Suppose further that 

(2.24) 

and that the bounds 
(2.25) 



< A 



A 



< \Dwj-i \ < AX in Qj 



hold for some A> 1. Then there exists a constant C4 depending only on n,p, ly, a, A 
such that 



(2.26) 



■j- \Du — Dwj \dxdt < C4 



^i-1 
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Proof. Let us begin by fixing several parameters appearing in the proof: 
1 . . „ , 1 



(2.27) 



7 



4(p-l)(n+l) ' ^ ■ ^ ■ " ■ 2(n+l) 

and, througliout the proof, we will apply Lemmas 12.21 and 12.31 with exponents 

Let Ci, C2, C3 be as in Lemma [2?2] and Lemma [231 respeetively, corresponding to 
these choices of a and e; therefore they ultimately depend only on nTp^v^a. We 
also set 



(2.28) 



A 



A 



In what follows constants denoted by c will only depend on v, a, A and will in 
general vary from line to line. We start to estimate the term on the left in (|2.26p 
with the aid of (|2.25p as follows: 



■f \Du- Dwj\dxdt < A(p-2)(i+7) i \Dwj_i\^P-^'>^^+-^^\Du^ Dwj\dxdt 

JQj JQj 

< c-f \Dwj - - Dwj\ dxdt 

(2.29) +C-/ \Dwj\^P~'^^'^^"^'''^\Du~ Dwj\dxdt. 

Appealing to Holder's inequality, together with (j2.19p and (|2.23p . gives us 
/ \Dwj - Dwj-i\'^P-'^^'^^-^''^\Du- Dwj\dxdt 

< ;^-(p-2)(l+7) ( I \Dw,-Dw,.i\^P-'^'^'+'>Uxdt 



(p-2)/(p-l) 



\Du - Dwi\P^^ dx dt 



< c(cr)4 



(l+7)(p-2) 



ciA 



1 



[(l+7)(p-2) + l](n+2)/[(p-l)n+p] 



[(l + 7)(p-2) + l](n + 2) ^ (p-l)(n + 2) 

(j> — l)n + p {p~l)n + p ~ 



But now, as 
(2.30) 

precisely for p > 2, (|2.24l) implies 

/ \Dw, - i5u;,_i|(P-2)(i+7)|^y _ dxdt < cM^'^ 

with c = c{n,p, a). Therefore (|2.29|) gives us 



\Du — DwA dxdt 



j 

(2.31) <c/ \Dw^\^v-'^)^^^'-i)\Du-Dw,\dxdt + c^-^^!^^^^ 
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We then continue to estimate the first term on the right in the above display. 
Applying Holder's inequality, together with (j2.2ip . and recalling that ^ = 1 + 27, 
we obtain for any a > that 

/ \Dwj\^P-^^'^^+'^'>\Du- Dwj\dxdt 
Jo 



< 



/ 

Jo 



,_^ {\Du\ + \Dw,\y-^\Du^Dw, 
{a+\u~Wj\)i 



1/2 



1/2 



i\Du\ + \Dwj \)P-^\Du- Dwj 1 2 
(a + \u — Wj\)^ 

\Dwj\^(P-'^\a+\u-Wj\)^dxdt 



dx dt 

dx dt 
. 1/2 



1/2 



< 



C2 
C-1 



1/2 



a 



(1-0/2 



1/2 



(2.32) 



1/2 



\DiVj\^'^P-'^\a + \u~Wj\)^ dxdt 



As the choice of a is still in our disposal, we set 
(2.33) 



Dwj\^'^P-^^\u-Wj\^dxdt 



for some small positive (5 G (0, 1) to get 



1/2 



\Dwj\^^P-'^\a+\u~Wj\)^ dxdt 



1/2 



< 2a«/2 ( 4 \Dwj\^^P~'^Uxdt 



2a«/2 



Notice that since Wj belongs to the parabolic Sobolev space LP{Tj \ VF^'^(Bj)) and 
^(p — 1) < p by (j2.27p . we have that a is finite by Holder's inequality. The presence 
in (j2.33|) of the parameter J, which shall be sent to zero at the end of the proof, 
guarantees that a is positive. In the above display, the integral on the right can be 
estimated by mean of Lemma 12.31 as 



/ 

JO 



DwS^P-^Uxdt 



<cf \Dwj-i - Dwjl^'^P-^Uxdt + cf \Dwj-i 

jQj Jo. 



^^P-^Uxdt 



1 MQj-i) 



Qj 

?(p-2)(ri+2)/[(p-l)ri+p] 



- cA«(f-2) < c , 
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owing to (|2.24p and (|2.25p . while the last constant c depends only on n,p,i',a,A. 
Thus (|2.32p . together with the last two displays, yields 

-I 1/2 



-/ \Dwj\^P-^^^^+^'^\Du~ Dwj\dxdt < c.f^ 
jQj V 



while applying Young's inequality together with an obvious estimation in turn gives 
Jq, 



(2.34) / \Du- DwJ\dxdt<^^- 

JQ, P 



for all /3 G (0, 1), where c = c{n,p, cr, A). Inserting this into (j2.3ip leads to 

again for all (3 G (0, 1), where c = c{n,p, v, cr, A) is in particular independent of j3. 
We then focus on a, which has been defined in (|2.33p . and split as follows: 

a < c 



(2.35) 



i^j \DiBj-i -Dwjf _^^.|?^2;<it^ 



u ~ Wj\^ dx dt 



S =: I1+I2 + S. 



By (ETUI), as C(p-l) = 9, we get by (I122)-(E1S1), together with Holder's inequality, 
that 

V (p-2)/9 / s 1/q 



h < c 



Wo_i — DwjV' dx dt 



\u — I'' dx dt 



1 



[(p-2)(n+2)+n+p]/[(p-l)«+p] 



Since 



'3-1 

{p-2){n + 2) + n+p 



2{p - 2) 



(p— l)n + p (p— l)n + p 

precisely for p > 2 we have, in view of (|2.24p . that 



> 1 , 



for c = c(n,p, u, cr). On the other hand, using condition (|2.25p we obtain 



u ~ Wj\^ dx dt 



Using Proposition 12. 11 estimate (|2.20p . and Young's inequality we get 

\ l/« / „ \ n/{n+l) 



\U — Wj 



1^ dxdt 



< 



u-Wj\^''+^'^''Uxdt 



< 



c -j- \Du — Dwj I dx dt sup / |m — wj \ dx 



n/(n+l) 



u — Dwj I dx dt 



^N-i 



n/{n+l) 



l/(n+l) 
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JQ 



\Du — Dwj \ dxdt + cr,_i 



ImKQj-ij 



Combining the estimates contained in the last three displays with (|2.35p leads to 



< c*rj_i -f \Du — Dwj \dx dt + c^,rj- 



^N-l 



with = c^,(n,p,i',a, A). Inserting this finally into (|2.34p with /? = l/(2c*) and 
then reabsorbing terms and sending 6 to zero finishes the proof. □ 

Next, a comparison estimate between Wj and Vj. 

Lemma 2.5. Let Wj and Vj be as in (j2.4p and (j2.5p . respectively, with j > 0. For 
every A > 1 there exists a constant C5 = C5{n,p^iy, L, A) such that the following 
holds: 

sup \Dwj \ + s < AX 

+ {\Dwj I + \Dvj\)P-^\Dwj - Dvj\^ dx dt 
HQ J 

(2.36) + / \Dwj - Dvj\Pdxdt < C5 [i^{rj)f . 

Proof. The result is based on the following estimate: 

j- {\Dwj\ + \Dvj\f-^\Dwj - Dvjf dxdt < c[uj{rj)]^ J- {\Dwj\ + sf dx dt 

that has been proved in [291 Lemma 4.3] and [J, with a constant c = c{n,p,v^L). 
The inequality in display p.36p follows by trivially estimating 

[^{'T])? i {\Dwj\ + sfdxdt < [Lo{r.j)f supd^w^l + sf < C5 [a;(rj)]^ 
HQ 3 hQ3 

and taking into account that p >2. □ 

Finally, along the lines of Lemma l2.41 we yet prove another comparison estimate, 
this time between u and Vj. 

Lemma 2.6. Let u be as in Theorcm \Ll\ and let Wj and Vj be as in (j2.4p and 

respectively, with j > 1. Suppose further that (j2.24p holds together with 



(2.37) 



sup \Dwj \ + s < AX 
iQj 

< \Dwj^i \ < AX in Qj 



for some A > 1. There exist positive constants cg = ce{n,p,iy, L, A) and cj 
CT{n,p,v,cr, A) such that the following inequality holds: 



(2.38) 



4- \Du — Dvj \dx dt < ceLd{rj)X + cy 
HQ, 



ImKQj-i) 

j.N-1 



Proof. We shall keep the notation introduced for Lemma 12.41 in particular, we 
recall (j2.28p . Inequality (|2.37D i allows to use (|2.36p so that 



/ ( I Dwj \ + \Dvj\)P-^\Dwj - Dvj l^dxdt 
HQ, 

(2.39) + / \Dwj - Dvj\P dx dt < cg [oj{rj)f . 
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With p' = p/{p — 1), by (|2.37p 9 we continue estimating as 

/ \Dwj - Dvj I dx dt < a(p-2)/p' I \Dwj-i\^P-^'>/P' \Dwj - Dv^ \ dx dt 



(2.40) 



|(p-2)/p' 



Dwj — Dvj I dx dt 



i I' 

HQ 3 



with c depending only on p and A and we arc using the notation in (|2.28p . As for 
the first term in the right hand side of (|2.40|) . we notice that since p > 2 we have 
(p — 2)/2 < [p — 2)/p' so that (j2.37P i allows to estimate 



/ \Dwj\'^P-'^'^/p'\Dwj - Dvj I dx dt 



< ^(p-2)V(2p);^(2-rt/2 i \Dw,\^P-^^/^\DWj ~ Dv,\dxdt 

so that, using Holder's inequality and (j2.39p yields 

/ \Dw^ I (P-2)/p' _ I dx dt 

< cA(2-p)/2 {\Dw,\ + \Dv,\)P-^\Dwj - Dvj\^ dxdt^ 

(2.41) < cA(2-P)/2c^/2u;(r,)Af/2 = ccj(rj)A 

for c = c{n,p, V, L, A). The second term in the right hand side of (j2.40p is estimated 
with the aid of Holder's inequality, (|2.23D . (|2.39p and finally Young's inequality, with 
conjugate exponents {p/2,p/ {p — 2)) when p > 2\ this means 



1/2 



/ \Dwj^i - \Dwj - Dvj I dx dt 



i/p' 



<c\4 \Dwj-i - DwAP^^ dxdt 



f \Dw 



i/p 



■.J - Dvj |P dx dt 



< a 



< c 



(p-2)/p' l/p 



1 



'j-i 



[(p-2)/p](p-l)(n+2)/[«(p-l)+p] 



1 



(p-l)(«+2)/[«(p-l)+p] 



A + uj{rj)X 



for c = c{n,p, V, L, a, A). Thanks to (|2.24p and (|2.30p . we get 

/ \Dw,^i~Dw,\'-P~^^/P'\Dw,-Dv,\dxdt<d Mgj^ 

with c = c{n,p,v^ L,(j^ A). By using the last inequality together with (j2.40p and 
(|2.4ip we conclude with 



+ uj(rj)\ 



■j- \Dwj — Dvj \dx dt < cuj{rj)X + c 
HQ 3 



for c = c{n,p,v, A) and c = c{n,p,i', L,a, A). Note next that the assumptions 
of this lemma fulfill also the assumptions of Lemma 12.41 Appealing then to (|2.26p 
and triangle inequality finishes the proof. □ 
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2.5. Comparison lemmas for SOLA. In this section we show that the basic 
comparison Lemmas 12.21 and 12.31 hold for SOLA in a suitably modified way. This 
fact ultimately allows to prove Theorem 11.51 by mean of the same proofs already 
given for the other theorems when considered for energy solutions. Before starting, 
since the setting here is the one defined in Theorem ll.51 we assume the existence of 
functions being local weak solutions to the equations considered in (|1.22l) and 
such that 



Uh 



and Duh 



Du a.e. 



(2.42) Duh e LP, Duh Du in LP~\ 

hold. The measures fih weakly* converge to /i and (|1.23p holds. Finally, since 
all the results we are interested in are local in nature, up to considering subsets 
compactly contained in Ht, we assume w.l.o.g. that the convergences in (|2.42p are 
valid in the whole fir- 

Moreover, in this section we keep the notation already introduced in Section [2. 21 
about the cylinders Qj , but the functions Wj will be defined in a different way: what 
matters here is not that they are solving Cauchy-Dirichlet problems as (|2.4p . but 
rather that they satisfy comparison estimates as in Lemmas 12.21 and 12.31 Indeed, 
once these lemmas are available for some suitably regular maps Wj then all the 
subsequent constructions can be replicated verbatim. More precisely, we have the 
following: 

Lemma 2.7. Let u be a SOLA to (jl.3|) as in Theorem ] L5\ and let e G (0, l/(n + l)]; 
moreover, let {Qj} be the sequence of cylinders considered in (j2.3p . Then there 
exists a sequence of functions Wj such that 

. e C'{T,;LUB,))nLUT,;W^^^{B,)) 

• Wj is a local weak solution to dtWj — div a(x,t, Dwj) ~ in Qj 



There exists a constant ci = ci(n,p, i/,e) such that 



(2.43) 




1/9 



Du - DwA'^dxdt 



< ciA 



IMOQ^-Jpar) 

A 



-| (n+2)/[(p-l)ri+p 



„N-1 



holds for any 0<q<p— l + l/{n + l)— e 



(2.44) 



• There exists a constant C2 = C2{n,p,v) 
{\Du 

Q, {a+\u-Wj\ 



f 

Jo 



DwAY-^\Du~Dw,\\ ai-« lA'KLQjJpar) 



{a ■ 



■ dxdt < C2 



C-1 



holds for any a > and ^ > 1 . 
• There exists a constant C3 = £3(71, p, e, a) such that the inequality 



(2.45) 



(2.46) 




1/9 



u — Wj I'' dx dt 



< csVj^iX 



l|/^|(LQj--lJpar) 

A 



{n+p)/[{p-l)n+p] 



holds whenever 0<q<p— 1+p/n — e and 

\ 1/9 



Dw 



'3-1 



Dwj \'^ dx dt 



< C3A 



1 lA^KLQ.-lJpar) 

A r^_i 



{n+2)/[{p-l)n+p] 



holds whenever 0<q<p— l + l/(7i+l)— s. 

Proof. The proof goes via approximation; we prove (|2.43p - (|2.44p . the proof of the 
remaining inequalities being similar. Let us fix a cylinder Qj and consider the 
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sequence appearing in (|2.42p and define 

be the unique solution to the Cauchy-Dirichlet problem 

{dtiifi — div a(x, i, Dwh ) = in Qj 
Wh = Uh on SparQj ■ 

Applying Lemma to this context gives 



(i 



1/9 



(2.47) 4 \Duh- Dwh\'^ dxdt] < ciA 



1 \^^h\{Q]) 



-| («+2)/[(p-l)n+p] 



for any 0<q<p— l + l/(n+l)— e and moreover 

(2.48 ■+ — — —7 dxdt<C2- 

Jq, {a+\uh-Wh\)^ C-1 

From (|2.47p and (|2.42p it follows that the sequence {Dwh} is bounded in U'~^{Qj). 
We now notice that Wh is an energy solution and solves an equation with Dini- 
continuous coefficients (sec Theorem 12.31 and subsequent Remark 12. ip . moreover all 
these equations satisfy assumptions (|1.4p - (|1.5p uniformly in h. Hence, by interior 
regularity theory (TB] and in particular by Corollary 12.11 it follows that the the 
maps {Dwh} are locally uniformly bounded in L°°. Again by interior regularity 
theory (see for instance the results in [52]) we have that the maps wt and {Dwh] 
are locally uniformly equicontinuous in Qj. Therefore, by Ascoli-Arzela's theorem 
and a standard diagonal argument, we may assume that, up to a not relabeled 
sequence, there exists a limit map w S J^focC^jS ^locfC^j)) ^^^h. that Wh ^ w 
in ifoc(^j' ^ioc'(-^i))' Dwh — > Dw and Wh w locally uniformly and almost 
everywhere. As a consequence, if weakly solves dtw — div a(a;, t, Dw) = in Qj . At 
this stage the proof of (|2.43p and (|2.44p follows taking wj := w, letting /i ^ cx) in 
(|2.47p - (|2.48p and using Fatou's lemma to deal with the left hand sides and p.23p to 
deal with the right hand sides. Inequalities in displays (|2.45p - (|2.46p can be similarly 
proved by approximation starting by the analogs of (|2.22l) - (|2.23l) . respectively, when 
written for Uh and Wh as already done for (j2.47p - (|2.48p . □ 



3. Proof of Theorems 11.11 and 11.21 

3.1. Proof of Theorem ll.il The proof goes in several steps and involves a rather 
delicate induction argument. In the following we select a Lebesgue point of the 
spatial gradient (a;o,to) € ^t, i-e., 

(3.1) lim + Du dx dt ^ Du{xo,to) . 

Almost every point in ilr, with respect to the Lebesgue measure in R"+-'^, satisfies 
such a property (see [131 Chapter 1, Page 8]). 

Step 1: Choice of constants and basic setup. With {xo,to) being fixed at 
the beginning we shall verify (jl.l9p with Q2r{xo,to) instead of Q^{xQ,to); this is of 
course causes no loss of generality. When finding constants c, Rq such that (|1.19p 
works, we choose positive numbers Hi, H2 appearing in the lower bound for A: 
(3.2) 

X>H,{4 i\Du\ + s)'^^Uxdt{^''''\H2 r \^m-o,to)) dg 



Q^A^o,to) J Jo e 
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where r G (0, Ro/2) and Rq is suitably small, to be determined in due course of the 
proof. The statement of Theorem II. II will be then proved with c := 2 max{i7i, H2}. 
Without loss of generality, we may assume that A is finite, since otherwise there 
is nothing to prove. In the rest of the proof certain constants will be deliberately 
chosen smaller/larger than necessary to emphasize the fact that their role is "to be 
very small/large". To begin with, having in mind to apply Theorems 12.31 and 12.41 
we set 

(3.3) A := 1000'*f^max{c2,C3,200}, B 10^ (5 10-^ e := 4'^^+^^ 

where C2 = C2{n.,p,i>, L) and C3 = 03(71, p, i>, L) are as in Corollarv 12.11 and The- 
orem [521 respectively. In particular, A depends only on n,p,v,L. Let <7i = 
ai{n,p,iy, L,uj{-)) and a2 = a2{n,p,v, L) be as in Theorems 12.31 and 12.41 respec- 
tively, both corresponding to the choices in (|3.3p . Set 



(3.4) c7:=min{c7i,a2,1000-P/^(p-l)-i/^A-(P-2)/Ar^^g-iVp| ^ (0,1/4). 

The choices made above guarantee that all fixed parameters A, B, e, 5, a depend only 
on n,p,i>, L,uj{-). Furthermore, let C4 and a be as in Theorem 12.51 corresponding 
to the choice of A in p.3p . In this way they both depend only on n,p, L. Next, 
let k be the smallest integer satisfying 

(3.5) C4cr'=" < — r and k>2. 

As all of C4, a and a depend only on n,p,v,L,uj{-), so does k = k{n,p,iy, L,uj(-)). 
We now proceed with the choice of Hi, H2 and Rq. We set 

(3.6) i?i := 100"/(P-i)lOV-4^. 

Then, taking ci = ci{n,p,v,ll {n + 1)) = ci{n,p,v) as in Lemma 12.21 C3 = 
C'i{n,p, V, L, 1/(71— 1), a) = C3{n,p, v,L,ijj(-)) from Lemma l^TSl and c^ = cr{n,p, v, L, 
U!{-), <t) = C7{n,p, V, L, from Lemma with A, a as in p.3[) - p.4p . we set 

(3.7) H2 := (2"cr-^('^+^)l06Amax{ci,C3,C7} 



Notice that since all the constants Hi and H2 arc built on ultimately depend on 
n,p,u,L,uj{-), we also have Hi,H2 = Hi, H2{n,p,v, L,uj{-)). We then pass to the 
choice of Rq. Looking then at Lemmas 12.51 and 12.61 and taking C5, cq corresponding 
to the choices of A and a just made in p.3p - p.4|) . we take R2 = R2{n,p, u, L, uj[-)) 
to be the largest positive number such that 

(3.8) ci/V-^('=+«) [u{R2)Y'^ + c-6a-(^^+i) £^ coig) f < ^ 

is satisfied. Finally, with Ri = Ri{n,p,v,L,ijj{-)) being as in Theorem 12.11 and 
Corollarv 12. 11 we let 

(3.9) i?o min{i?i,i?2}/4. 

From now on, with A as in (|3.2p . 2r < Rq and (32^(^07^0) C fix, we adopt and 



fix the basic set-up described at the beginning of Section [Z2l We therefore denote 
rj = a^r for j > 0, Qj = Qrj{xQ,to) and the comparison maps Wj and Vj defined 
in (|2.4p and (|2.5p . respectively. 

We finally record a few immediate consequences of the choices done so far for 
Hi, H2, Rq. Observe that in (j3.2p the latter potential term may be further estimated 
from below as follows (recall that rg — r): 

f^- \fi\{Ql{xo,to)) dg 

Jo 9"-' Q 
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> 



E 



fi\{Q^^{xo,to)) dg , f^^« |M|(g^(a-o,to)) dg 



,JV-1 



(3.10) 



E 

1=0 



+1^ 



dg , ImKQo) 



(2ro 



log 



E 



log 2 
9W-1 



ro e 

>l(Qo) 



SO that we conclude with 
(3.11) A > i/i2-^/(P 

The choice of Hi in 
(3.12) 



" (/ 'I 



Du\ + s)P-Uxdt 



„JV-1 



i/(p-i) 



> 



-"E 



guarantees for instance that 
A 



s < 



600 



Notice that in the decomposition p.lO[) we can always assume that 

(3.13) |//|(aparQO = 

for every i > 0, since this amounts to change the integrand of the first integral in 
(|3.10|) in countably many points. The choice of H2 in (j3.7[) together with (j3.13p . 
allows to deduce that 



E 



E 

i=0 



< 



< 



l/^KLQdpar) 



^JV(fc+4) 



2^10Mmax{ci,C3,C7} 



[n(p-l)+p]/(n+2) 



-A < A. 



^^'^^^ " 2^106Amax{ci,C3,C7}' 

Notice that in the second estimate above we have used the fact that, since p > 2, 
then n{p — 1) + p > n + 2, and that the quantity in brackets is smaller than 1. 
Recall now the choice of Rq in (|3.9|) and observe that, if r = tq G (0, i?2/4], then 



R2 



K.)^ - Ef -(^)- 



dg 

uj{g} — 
9 



(3.15) > logf-)^c^(r,+i)+log4L.(ro) >a^c^(rO 

holds, so that p.8p implies 



j=0 



i=0 



(3.16) 



_l/p 7V(fe+4)r 



2/p _ 



C6 



< 



2^106 



Remark 3.1. When dealing with equations of the type in p.20p . i.e., equations 
without dependence on x, Theorem 12.31 is not any longer needed and therefore 
no dependence on a;(-) occurs in the constants. Actually, we do not need the 
comparison functions Wj and we can just use Vj. 

Remark 3.2. An ambiguity occurs when considering (j3.13p . since p.l3p is auto- 
matically satisfied when ^ S L^, which is precisely the case when proving Theorem 
11.11 Here we nevertheless wanted to emphasize the fact that we can assume (|3.13p 
at this stage also when ^ is genuinely a measure. This fact will play a role when 
proving Theorem 11.51 below. 



RIESZ POTENTIALS AND NONLINEAR PARABOLIC EQUATIONS 



27 



Step 2: Exit time argument. After having fixed the relevant constants in the 
previous step, we consider, for indexes i > I, the quantities 



/ , \ i/(p-i) 



l^j \Du\P-^dxdtj +2a-^ -j- \Du~ {Du)Q^\dxdt. 

By (inH) and ([XTT]) it follows that 

M f r \i/(p-i) 
Ci<&a-^-—^ [J- \Du\P-^dxdtJ < 



1000 



Furthermore, without loss of generality, we may assume that there is an exit time 
index Ze > 1 such that Ci > A/1000 whenever i > ie and Ci^ < A/1000. Indeed, if 
such an index does not exist, we have - in view of Ci < A/1000 - a subsequence 
of indexes such that d- < A/1000 as j — >• oo. But, as we assume that {xo,to) 
is a Lcbesgue point of Du, then also 



(3.17) \Du{xo,to)\ = lim \iDu)Q^,\ < limsup C,, < 



A 



1000 



would hold and the proof would be complete. Thus, from now on, we shall work 
under the assumption 

(3.18) >^ for *e{*e + l,V + 2,...}, and ^ < ■ 

Note that in (j3.17p we are using a limit computed on a sequence of shrinking 
intrinsic cylinders, therefore different from those considered in (|3.1|) . On the other 
hand, define the set C\ as 

Cx := < {xo,to) S fir '■ liui -f Dudxdt = Du{xo,to) > 

for A > 0. Basic properties of maximal operators - see for instance [43l Chapter 
1, Page 8] - imply that this set is actually independent of A and, in particular, 
C\ = Ci =: C for all < A < c« and therefore p.l7p is completely justified in view 
of the assumed property in p.l[) . 

Step 3: Induction scheme. In order to prove (|1.19l) using p.lSp . we apply an 
induction argument. To shorten the notation, we set, for i > 0, 

(3.19) E,:=J^ \Du- {Du)Q^\dxdt, a, := \{Du)q^\ . 

We shall consider, in our iterative setting, the following conditions: 

/ \ / \ i/(p-i) 

Indi(j) : if \Du\P-Uxdt] + i -L \Du\P-^dxdt] <A 

\JQj-i J \Jqj J 

for given integer j > and 

Ind2(.): t i^.<^E^« + i^E-(-)^+5 £ 

i— ie+l i—i^ i—it^ — 1 ^ 

for given integer j > if.. The constants cg, cj are those defined in Lemma 12.61 with 
the choices of A and a made in p.3p and p.4p . respectively; in this way, they 
ultimately depend on v, L, cj(-). 
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The goal is now to show, by induction, that Indi(j) holds for every integer 
j > ie and Ind2(j) holds for every integer j > ig. Indeed, this will immediately 
prove Theorem 11.11 as 



\Du{xo,tQ)\ = lim Gj < limsup 4 jDuj^ ^ dxdt 



i/(p-i) 



< A. 



Now let us remark that Indi(ze) is automatically satisfied since 
/ „ \ 



(3.20) 



\Du\P-^dxdt 



A 



1000 



Therefore, the rest of the proof develops according to the following scheme: 



(3.21) 

(3.22) 

and 

(3.23) 



Indi(ie) 



f Indi(j) 
1 Ind2(j) 

Indi(j) 
Ind2(j + 1) 



> Ind2(ie + 1), 
Ind2(j + 1) yj>ie 

Indi(j + 1) \/j>ie 



We remark that in the following, unless otherwise stated, whenever we are con- 
sidering Indi(j) we will do it in the general case j > ie, while, when considering 
Ind2(j) we will do it for all indexes j > ie- 

Step 4: Upper bounds implied by Indi (j). Towards the proof of the induction 
step, here we assume that Indi(j) holds for a certain index j > and exploit a 
few consequences of this. With the integer k being defined as in p.Sp . observe that 
Lemma and p.l4p imply that whenever I S {0,l,...,fc-|-l} the following holds: 



\Du- Dwi-i\P~'^dxdt 



Qj-i+i 



< 



i/(p-i) 



Du - Dwj-i\P~^ dx dt 



(3.24) 



1 WQi-i) 



-| (ri+2)/[(p-l)ri+p] 



< 



2^106 



A. 



Similarly, the inequalities 
(3.25) 



i/(p-i) 



Du - Dwj\P^^ dxdt 



- 2^106 



hold as well for Z G {0, 1, . . . , fc + 1}. Using with I = and Indi (j) we get 

\ i/(p-i) 
\Dwj.i\P-^dxdt] 



f 



< 



2^10' 



:A + 



l/(p-l) 



dx dt 



< 2A. 
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By Corollary 12. II (applied to wj-i) we then obtain, also recalling p.l2p . that 

(3.26) sup iDwj.il + s < C2(A + s) + C2A2-p -/ {\Dwj^i\ + s)p^^ dx dt < A\ . 

In a completely similar way, by using this time (|3.25p instead of p.24p . and again 
Indi(j'), we can prove 

(3.27) [j \Dwj\P-^ dxdt\ < 2A 
and then 

(3.28) su^\Dwj\+ s < A\. 

Furthermore, keeping (j3.16p in mind, estimate (|3.28p and Lemma 12.51 allow us to 
deduce 

i/(p-i) 



(3.29) \Dw.j-Dvj\P-^dxdt\ < cg/^ [w (rj )]'/'' A < 



Ar(fe+6) 



2^106 

so that, for Z 6 {1, . . . , fc + 1} it holds that 



y' Qj + l 



\Dwj ~ Dvj\P-^ dxdt] < ,,,_^ A. 



2^106' 



Combining the above estimate with p.25p gives 



i/(f-i) ^ 



(3.30) if \Du - Dvj\P-^ dx dt 

for alU G {1, . . . , fc + 1}. Next, observe that (|3^ and ((3?29l) imply 

\ i/(p-i) 



i 



\Dvj\P-^dxdt 



/ \ / X i/(p-i) 

< -/ \Dwj \P-^ dxdt] + f \Dwj - Dvj \P-^ dx dt ) 

\JhQ. J VhQ, J 

<2^+iA + A<4^A 

so that Theorem O and ((XT^ finally yield 

(3.31) snp\Dvj\+s<C3{X + s) + C3X'^-P + {\Dvj\ + s)p-'^ dx dt < AX . 

Step 5: Lower bounds implied by Indi(j). Here we still exploit a few con- 
sequences of assuming Indi(j) for some j > i^. In particular, we derive suitable 
lower bounds for Dwj-i and Dvj. First, let us first show a few oscillation reduction 
estimates. Since we have already established the upper bound for Dwj-i in p.26p . 
Theorem l2 . 31 applied to Wj-i (with S = 10^^ as in (j3.3p and recalling also the choice 
made in p.4p that guarantees cr < ci so that Qj C ctiQ^-i) gives 

(3.32) oscDwj^i<^. 

For Dvj we instead have p.3ip and hence Theorem 12.51 and p.Sp imply 

osc Dvj < CiCT°'^X < —A. 
Q,+. ^ - - lOe 
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Using this together with (|2.2p and p.30p gives 
2(7"^/ \Du- {Du)Q^^Jdxdt 

<4a"^4 \Dvj - {Dvj)Q^^Jdxdt + Aa-^ -I \Du-Dvj\dxdt 

^ Qj+k ^ Qj+k 



< 4cr~^ osc Dvj + 4ct"^ / \Du- Dvj \dxdt < ^ 



10^ 



But, as Cj+A; > A/1000 for j > Zg, we then necessarily have that 
(3.33) J \Du\P~Uxdt 1 > 



' Qi + m + fc 



2000 



Next, again by using triangle inequahty and (j3.24[) (for Z = fc, fc + 1), we have 

t \Du\P-Uxdt\ 

/ X i/(p-i) 

(3.34) ^T7i^+ E f |i?^«,-iridxdt 

so that, as fc > 2 and p.33p holds, we also get 

/ \ i/(p-i) 

2sup\Dwj-i\ > V f \Dwj-i\P-'^dxdt] 

^^■^^^ - 2000 106 - 104 ■ 



Arguing as for p.34p . and using p.30p . we this time have 

i-'i'-' 



/ , \ 
_ \Du\P~^dxdt\ 

= _1 V-^Qj + m + fc / 

/ \ i/(p-i) 

< ^+ Y if \Dv,\P-^dxdt 

m=-l + 

SO that, using again p.33p . and recalling that fc > 2, we conclude with 

i/(p-i) 

^ 10* 



(3.36) 2sup|i:)wj|> ^ ( -/■ dxdt] > ^'^ 



The inequality in display p.35p yields the existence of a point (i, t) G such 
that iDwj-i (a;, t)| > A/IO** and therefore the oscillation control in p.32p . recalling 
also p.26p and that Qj C ■^Qj-i, gives 

(3.37) ^<JL<\Dwj.^\<A\ mQ,. 
Finally, (|3.36p and p.3ip . and the fact that Qj+i C jQj, give 

(3.38) — = ^ < sup \Dvj \ < sup \Dvj\ + s<A\. 
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Step 6: Further consequences of Indi(j). We again assume Indi(j) for an 
arbitrary chosen j > i^. Now, on one hand p.l4p . p.28p and p.37p ahow to apply 
Lemma 12.61 and obtain 



(3.39) 



■f \Du — Dvj\dx dt < CQU){rj)\ + 



On the other hand, recalhng that aQj = Qj+i, Theorem 12.41 with e = 4 (^+4) is 
at our disposal by (|3.38p and it yields 

(3.40) ■/ \Dv,-{Dv,)Q^^,\dxdt<i'^^+^U^ \Dv,-{Dv,).Q^\dxdt. 

Now, by (|3.39p . recalling the definitions in p.l9p and using (|2.2p repeatedly, both 
4 \Dvj - {Dvj)iQ^\dxdt 

< 22^+1-/ \Du- {Du)Q^\dxdt + 2^+^ -f \Du~Dvj\dxdt 

Jq, ' JhQ. 



< 2^^+'E, + 2^+2c6w(r,)A + 2^+^ 



and 



2-1 \Dvj~{Dvj)Q^^,\dxdt 
JQj+i 

> -I \Du-{Du)Q^^,\dxdt-a-^ -I \Du - Dv, 

ImKQj-i) 



dx dt 



> Ej j^i — G c6w(»'j)A — (7 c^ 



^3-1 



hold. Combining the inequalities in the last three displays gives 

l/^KQ.-i)' 



(3.41) 



1 2ce , , 2cr 



j,N-l 



Step 7: Verification of Ind2(ie + 1) and Ind2(j' + 1). Here we prove (|3.2ip 
and p.22p . The outcome of Step 6 is that p.4ip holds whenever Indi(j) holds, for 
every j > ie, therefore, since Indi(ie) holds as established in (j3.20p . taking j = ig in 
(I3.4ip yields Ind2(je + !)• Wc now come to the proof of (j3.22p . that is the validity 
of Ind2(j + 1); this simply follows by summing (|3.4ip . that holds by the assumed 
validity of Indi(j), to the inequality yielded by the definition of Ind2(j), which is 
also assumed in p.22p . 

Step 8: Bounds for aj and Ej. It remains to prove p.23p . For this we now 
assume that Indi(j) and Ind2(j + 1) hold for some j > ie and derive bounds for 
the quantities defined in p.l9p . By Ind2(j + 1) and easy manipulations, we note 
that (EH), and < A/1000 imply 

}_^E, < 2E,^ + -^}_^uj{n)X + -^}_^—^ 

i=i^ i=0 i=0 t 

a^X . ^"^A . G^X 



(3.42) 



< 2K, 



1000 



<g'^C, 



< 



1000 - 500 
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Using this together with the obvious estimation (vahd whenever i > 0) 
Oi+i — a; < -+■ \Du — {Du)q.\ dx dt 



< 



^/ \Du-{Du)Q^\dxdt^a-^E, 

Jo, 



\Qi+l\ JQ, 

we get, after telescoping the previous inequahties and using (j3.42|) . that 

^ 1000 ' 500 - 200 



(3.43) «m<a.+<^-"E^^=^T4l + 4^ ' 



1=1, 



for aU I G {ie, ■ ■ ■ Here we also used that, thanks to p.lSp . we have 
(3.44) < [ 4 \Du\P-^dxdt] <G,^ < 



1000 



step 9: Verification of Indi(j + 1). Here we prove p.23p . thereby concluding 
the proof. We actually have to prove that 

/ \ i/(p-i) 

(3.45) T \Du\P-Uxdt\ <X. 

J71 = -l ^Qj + m + l j 

To this end, we estimate using (j3.24l) (with Z = 1, 2) as follows: 

\ i/(p-i) 
\Du\^-^dxdt\ 



til 

■n=-l V-^Qi + ^ + l 



/ \ i/(p-i) 

< Y if \Du- Dwj^i\P-^ dxdt\ 

m=-l V-^Qj + m + i / 

/ \ i/(p-i) 

E 4 \Dwj^i\P-'^dxdt\ 

m=-l \JQ] + m+l ) 



(3.46) ^Wm^^^ (/ \Dw,^irUxdt 

71— 1 VQj + r^^+l 



7n—- 

/ , \ 



We further estimate the latter term on the right hand side by simply applying 
triangle inequality as follows: 



E 4 \Dwj_i\P-^ dxdt 



< f \Dw,.,r'\{Du)Q^^„^^Adxdt 



+ f \Dwj-i\P-^ {\Du^ {Du)q^^^^,\ + \Du- Dwj_i\) dxdt. 

Using Young's inequality with conjugate exponents {p — l)/{p — 2) and p — 1 (only 
when p > 2) we get 



E f \Dw,^^r'\{Du)Q^^^^,\dxdt 

^ — 1 J Q■i-\-m .-\-^ 



m— — 1 
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<^ V 4 \Dwj-i\P-'^dxdt 

^ m=-l 

Matching the inequahties in the last two displays, reabsorbing terms, and us- 
ing dXMl), yields 

» 

m=-l "^Qj + m+i ?ri=-l 


+ l)(AAf-2 ^ J. (|DM-(Dw)Q^^,„^J + |i:)ii-L>u;j_i|) dxdt. 

Notice that we have used that Qj C ^Qj-i in order to apply p.26p . Now we 
estimate all the terms in the right hand side of the above inequality. Using (|3.43p 
and p.44p (this last one only when j = ie) from Step 7 we have 

m— — 1 ^ ^ 

while using p.42p gives 

f \Du-{Du)Q^^^^,\dxdt = E,+,+E,<—-. 

Finally, using (|3.24p (with I — 1,2) and Holder's inequality yields 



E 7 \Du-Dwj^i\dxdt<-^—^\. 

Connecting the inequalities in the last four displays, and recalling the very defini- 
tions of A and a in p.3p and p.4p respectively, gives us 

^ 4 \Dwj-i\P-^ dxdt < ' 

Inserting the last inequality into p.46p leads to p.45p . Therefore p.23p is verified 
and the proof of Theorem 1 1.11 is complete. 

3.2. Proof of Theorem 11.21 With the standard cylinder Qr = Qr{xo,to) C Ht 
being fixed in the statement, let us consider the function h(X) := A — cA(X) where 



i/(p-i) 



P-2 



A{X) := X— 



^— [ {\Du\ + s + lY-^dxdt\ +l>'^^[xoM:r) 

i-r\ Jq>^ I 



:= 1^-^ {\Du\ + s + l)P-Uxdtj +l5'^(xo,to;0 

and c > 1 is again the constant appearing in Theorem 11.11 it depends only on 
n,p, V, L, aj(-). Here it is = Q^{xo,to). We are actually considering the function 
h{-) to be defined for all those positive A which arc such that C ^It', observe that 
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since p > 2 then the domain of definition of h{-) includes [1, cxd) as C Qr C ^It 
when A > 1. Moreover, observe that again when A > 1 we have 

(3.47) Ai\)<\^i^j {\Du\ + s + l)P-Uxdtj +I'^{xa,to;r). 

The function h{-) is obviously continuous and moreover h(l) < since c > 1 and 
^(A) > 1. On the other hand, p.47|) implies that h{X) — J- oo as A — >■ oo. It follows 
that there exists a number A > 1 such that h{X) = and therefore A = cA{X). In 
particular, A satisfies (|1.19p . Therefore we can apply Theorem 11.11 that gives 

A + \Du{xo, to)\ < 2A = 2c^(A) . 

Using in turn Young's inequality with conjugate exponents (p — 1, {p — l)/(p — 2)) 
when p > 2 and p.47|) . we have 

2cA(A) <^ + c-f {\Du\ + s + 1)P-^ dx dt + 2cI^{xo,to;r) 

where c depends only on n,p,iy, L,uj{-). The proof follows connecting the inequali- 
ties in the last two displays. 

4. Proof of Theorems 11.31 and 11.41 

4.1. Proof of Theorem 11.41 The proof consists of several steps; some of the 
arguments of the proof of Theorem 1 1 . 1 1 will be re-proposed in order to, this time, 
control the degeneracy rate of the equation. 

Step 1: Basic setup and smallness conditions. Since we are assuming that 
I^(xo, to; r) is locally bounded for some r > then by Theorcm ll.2l we have that Du 
is locally bounded in fi^, too. Moreover, since we are proving a local statement, up 
to passing to open subsets compactly contained in ilx, we can assume w.l.o.g. that 
the gradient is globally bounded, therefore letting 

(4.1) A := ||Z?u|Ii^(n^) +S + K «). 

From now on our analysis will proceed on cylinders of the type C ^It- We shall 
prove that for every s £ (0, 1), there exists a radius = r^{n,p, v, L, a;(-), /i(-)i s) > 
such that 

(4.2) E{Du,Q^)^-f \Du-{Du)Qx\dxdt<\e 

holds whenever g S {0,r^] and Qg C Ht- Once this fact is proved the VMO- 
regularity of Du in the sense of Theorem 11.41 follows by an easy change-of- variables 
argument as A is now fixed in (|4.ip . 

Towards the application of Theorems I2.3ll^^ and with C2 = C2(n,p,i>, L) and 
C3 = C3{n,p,iy,L) being as in Corollary [^TT] and Theorem l2.21 respectively, we start 
fixing 

10004p^max{c2,C3,200} 10^ e _ s 

4.3 A := , i> := , d := — - , e := — rirr . 

^ ' e £ ' 105 ' 10^^ 

With the choice in (|4.3p we determine ci = ci (n,p, i^, i, a;(-), e) and CT2 = f2(n,p, 
iy,L,e) from Theorems 12.31 and 12.41 respectively. Next, we fix 

(4.4) cr := min{cri, era, 16"^} G (0, 1/4) 

and notice that a only depends on n,p,iy,L,uj{-)^e. Such choices, looking at Lem- 
mas [221 [23] and determine the constants ci, C5, cq and C7 as depending only on 
n,p, v, L, uj{-), e. We again need some limitations on the size of the radii considered; 
here i?i = Ri{n,p, v, L, still denotes the radius considered in Theorem 12. II and 
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Corollary 12. II Moreover, we select a new radius i?3 = i?3(?i,p, z^, L, cj(-), /^(•), e) in 
such a way that the following smallness conditions hold: 

(r,+2)/[(p-l)n+p] 



(4.5) 

and 

(4.6) 



sup sup 

o<e<R3 {x,t)GQ.T 



,N-l 



Xg 



< 



lO^CiCsCeCy 



and we this time set Rq :— min{_Ri, i?3}/4; as a consequence we have Rq = 
Ro{n,p, V, L, uj{-), fi{-),e). Notice only that the possibility of gaining (j4.5p follows di- 
rectly by the assumption (jl.21|) . and from the fact that, since A > 1 and p > 2, then 
Q^{x,t) C Qg{x,t). We now fix a cylinder Q = Q^{xo, to) C fir with r G {uRq, Rq], 



and accordingly to the setup defined in Section 
parabolic intrinsic cylinders as follows: 

(4.7) Q, = Q^(xo,to) ^ 



define the chain of shrinking 



(T-'r, where rG{aRo,Ro] 



for every integer j > 0. The related comparison solutions Wj and Vj are accordingly 
defined as in (|2.4p and (|2.5p . respectively; finally, we denote 



j 



\Du — {Du)q. \ dxdt . 



We shall preliminary prove that 

(4.8) Ej+i < Xe Vj e N n [1, oo) . 

Now, it is obvious that if for a given j > 1 it occurs that 



(4.9) 




Du\P-^dxdt 



< 



Xe 
50 ' 



then (|4.8p holds, therefore we can confine ourselves to assume that (|4.9p does not 
hold. To prove (|4.8p in this last case we start proving the following implication, 
which is valid whenever j > 1: 



i/(p-i) 



\Du\P-^dxdt 



> 



' Qj+i 



(4.10) 



E, 



< 



300 



:E, 



Ae 
50 

2C6 



Indeed, as in any case we have Ej < 2X and (j4.5p - (j4.6p hold, it follows Sj+i < eA/50 
as a direct consequence of (j4.10p when (|4.9p does not hold. Notice that wc have 
used that (n + 2)/[(p — l)n + p] < 1, which holds since we are considering the case 
p > 2. Therefore to prove (|4.8p we are reduced to check the validity of (|4.10p . 



Step 2: Proof of (|4.10p . Wc note that the setting of the shrinking cylinders {Qj} 
and related comparison maps Vj , Wj defined in (j4.7p is completely similar to the one 
adopted in the proof of Theorem 11.11 only the choice of the constants differs. To 
prove (j4.10p is in turn sufficient to prove the following group of inequalities for the 
choices made in (|4.3p : 



(4.11) 



sup \Dwj \ + s < AX 



(4.12) 



4 ^ 

A ■' 



-i| < \Dw,^i \ + s < AX 



in Qj 
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A 

< 



(4.13) — < sup \Dvj \ < sup \Dvj \ + s < AX . 

13 - 



Indeed, taking these for granted, let us sec how to conclude with the proof of (|4.10p . 
Inequalities in displays (j4.1ip - (j4.13p arc completely analogous to (j3.28|) and (j3.37|) - 
p.38p . respectively. We can therefore exactly argue as in Step 6 of the proof of 
Theorem O using (|ilT|) - (|il^ we get ([535)1 by Lemma while allows 
to use Theorem 12.41 with e = lO^^^e thereby yielding 



/ \Dv-i - {Dvj)q^^, \dxdt< f \Dvj - {Dvj)iQ^ I dx 



dt 



which plays the role of p.40p in this context. Proceeding as in Step 6 we finally 
arrive at the inequality in display (|4.10p . which is the analog of p.4ip in this context. 
It remains to prove (j4.1ip - ()4.13p . Using Lemma and (j4.5p . we have 



\ i/(p-i) 

\Du - Dw j-i\P~^ dxdt I 



(4.14) <cr"^(/ \Du~ Dwj^i\P'^ dxdt] < 



i/(p-i) 

^ 104 



i 



By (|4.14p . keeping (|4.ip in mind, and using triangle inequality we get 

\ i/(p-i) 
\Dwj^i\P-Uxdt] <2A 

so that, applying Corollarv l2.1l to Wj^i wc finally come to 
(4.15) sup + s < sup \Dwj-i\ + s < AX 

that proves the right hand side inequality in (j4.12p . In a completely similar way 
we obtain also (|4.1ip . We now prove the left hand side inequality in (|4.12p . Again 
using (|4.14p . we get 

sup|L»u;j_i| > 4 \Du\P~^dxdt\ 



I ^ ^ I ^ 1 . . \ Xs Xc Xs 

\Du-D^,^,r'dxdt\ >^-T5I>n5o- 



The lower bound in the last display provides the existence of a point {x,t) S Qj 
such that \Dwj^i{x,i)\ > Ae/200 while (|4?T5)) allows to apply Theorem [231 (to 
Wj-i in iQj-i) getting, thanks to the choice of 5 in (j4.3p . 



Xe_ 



o^c Dwj^i < 



The last two inequalities finally give \Dwj-i\ > Ae/10^ in Qj. Summarizing, re- 
calling the choice of A in (j4.3p and again (|4.15p . the proof (j4.12p turns out to be 
complete. It remains to prove (|4.13p . Using (|4.1ip . Lemma [^31 and finally (|4.6p we 
have 



V 



\Dwj - Dvj dxdt] < [w(rj)]^^^ A < 



106 
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Combining the above estimate with (|2.19p and (|4.5p gives, after a few standard 
manipulations, 



(4.16) 





Dvj \P-^ dxdt] A 4 \Du^ Dv, dx dt 



1 WQj) 



-| (ri+2)/[(p-l)ri+p] 



A r 

3 



N- 



< 



106 - 105 



Therefore, using (|4.16p . Minkowski's inequahty and recalhng the definition of A in 
(|4.ip we have 

\ i/(p-i) 




DvAP-'^dxdt] <A + ^<2A 

J I I - 10^- 



so that the right hand side inequahty in (|4.13|) follows by Theorem 12.21 As for the 
left hand side, notice that the first inequality in (|4.10p and (|4.16p imply 

sup \Dvj \ > ( -f \Du\P-^dxdt ' 




Du-Dv.F-Uxdt\ >^o~^^Wb 



"J 

'Qj+i 



so that (|4.13p is completely proved. 

Step 3: Interpolation of radii. With (|4.8p at our disposal we can finally con- 
clude the proof of Theorem 11.41 by letting := o^Rq. Indeed, consider g < a^Ro; 
this means there exists an integer m > 2 such that cr™+^^o < P < (j'^Rq. There- 
fore we have g = cr^r for some r G {<jRo, Rq] and (|4.2p follows from (|4.8p with this 
particular choice of r. 

4.2. Proof of Theorem 11.31 The proof is now based on a combination of the 
arguments of Theorem 11.41 with those which are more typical of the elliptic case; 
we report everything in full detail for the sake of completeness and readability, and 
also because a certain number of modifications is really needed. We shall therefore 
keep the notation introduced in Step 1 of the proof of Theorem II. 41 Essentially, we 
are going to use the same choices in (|4.3p - (|4.6p but using an additional smallness 
condition on the radii used; in this way we can use both the inequalities in the proof 
of Theorem [131 and the result of Theorem ll.4l We consider a cylinder Qq g fix and 
prove that for every e > there exists a radius < c?par(Qo7 5r2T)/2, depending 
only on v, L, a;(-), fJ.{-), e, such that 

(4.17) l(^w)Q^(:i.o,to) " iD^)Q^{xoJo}\ ^ liolds for every g, p e {0,r^] 

whenever [xq, to) G Qo- This proves that Du is the local uniform limit of continuous 
maps - defined via the averages - and hence it is continuous. The rest of the proof 
goes in two steps. 

Step 1: Dyadic sequences and continuity. To begin with the proof of (|4.17p 
we recall that Ri is determined in Theorem l2.1l Moreover R3 is determined in (j4.5p - 
(j4.6p with the constant ci , C5 , cg and cj obtained Theorem 11.41 and corresponding 
to the choices made in (|4.3p and (|4.4p . Next, we take yet another positive radius 
R4 < c'par(Qo7 5fiT)/2 such that 

(4.18) sup / -j^ + / uj{g) — < 

{x.t)eQo Jo ^" 



Ag^-i g Jo ' ' Q - IO6C6C7 
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and, recalling the definition in (|2.16p 



(4.19) sup sup E{Du,Q){x,t)) < . 

Q<Q<Ri {x,t)eQa 

Observe that it is possible to make the choice in (j4.19p thanks to Theorem 11.41 
and in particular to (j4.2p . that ensures that i?4 can be chosen in a way that 
makes it depending only on n,p, i/, L, w(-), //(•), £. Finally, this time we set i?o '■= 
min{i?i, i?3, i?4}/2 and take everywhere r < Rq so that i?o ultimately depends 
again on n,p,iy, L,lo{-), fi{-), e only. The sequence of shrinking cylinders {Qj} is 
now defined as 

Qj = Qr, {xo,to) , rj = (T^Ro, for j > , 

while A > 1 is still defined as in (|4.ip and a in (|4.4p . By (j4.18p . computations 
similar to those in p.lOp and p.lSp then give 



=0 



In Step 2 we will prove that 

(4.21) \{Du)q^ - {Du)q^ I < — holds whenever 2<k <h. 

Here we show how to use (|4.2ip to finish the proof and to verify (|4.17p with the 
choice := ct^Rq. Indeed, let us fix < p < p < r^. This means that there exist 
two integers, 2 < k < h, such that a'^+^Ro < g < a^Ro and a'^+^Ro < p < a^Ro- 
Applying (I4.19P we get 



\i.Du)QX( to)-{Du)Q^^^\ < 



i \Du- {Du)Qx^^^^t^-j\dxdt 

Qk+l 



< 



— . — f \Du- [Du)QX,^^^^t^)\dxdt 

\Qk+l\ JQ^Ixo,to} 



and; similarly. 



< a~'^E{Du,Q^^{x„,to)) < 



Xe 
10 



Using the inequalities in the last two displays together with (|4.2ip and triangle 
inequality establishes (|4.17p and the proof is complete, modulo the content of the 
next and final step. 

Step 2: Proof of (|4.2ip . The preliminary observation to make is that, with the 
choices made here, (j4.10p holds in this setting as this ultimately relies on (j4.5l) - (|4.6p . 
that indeed are in force by the choice of Rq. To continue, let us consider the set C 
defined by 

C:^ {j eN : \ 4 \Du\P-'^ dx dt] < 




50 

and, accordingly, we then define the sets 

Cr = {j eN ■.i<j<i + m, i(zC, j (^C\ii> i} 

for TO G N and, finally, the number je := min C. Note that it may happen that 
je = oo; in this case C is empty and the first inequality in (j4.10p holds for every 
j > 1. The idea is now to employ (|4.10p on suitable sets C™ using the indexes i as 
a sort of exit time indexes; the difference is that they can be countably many now. 
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We can now prove (|4.2ip . obviously assuming k < h. The first case we analyze is 
when k < h < je', we use (|4.10p and the definition of je to infer that the inequality 



(4.22) E,+,<lE, + ^u{r,)X + ^ 



holds for every j G {fc — — 2}. Summing up the previous inequalities easily 

yields 

j:=fc j=0 j=0 j 

where we have used (|4.19l) - (|4.20p . therefore (|4.2ip follows since 

h-l 

\{Du)q, - iDu)Qj < J2\^^''h.+i - iDu)Q^ 

i=k 

h-l 

< V4 \Du- {Du)Q.\dxdt 

^ \Du-{Du)Q^\dxdt 

(4.23) = .-"E^^^^^- 

i—k 

The second case we consider is when je < k < h, where we prove (|4.2ip through 
the inequalities 

(4.24) \(^Du)qJ<^ and \{Du)qJ < ^ . 

In (|4.24p . we prove the former, the argument for the latter being the same when 
k > je, otherwise \{Du)qi,\ < Ae/25 is trivial if A: je G C. If h £ C, the first 
inequality in (|4.24p follows immediately from the definition of C. On the other 
hand, ii h ^ C, then, as ft. > je, it is possible to consider a set C™*" with rrih > 0, 
such that h e C""; notice that h > a.s h ^ C 3 ih- Then (jil^ gives that 
holds whenever j G {ih, . . . ,ih + mh — 1}. Summing up and performing elementary 
manipulations gives 

ih+nih ._oo A- °° \ \ 2N \ 

\^ F<9P I ^ , . , 4c7 ^ Ia^KQ,) ^ g'^Ae 
2^ ^, < 2S,, + ^ 2^ w(rj) A + ^ 2^ w-i < 

i=ih j=0 j=0 j 

where again we have used (|4.19p - (|4.20p . Therefore, as in (|4.23p . we have 

R, < 

50 



and then, using that \{Du)q.^ \ < Ae/50 as G £, we have 

\{Du)q,\ < \{Du)qJ + \{Du)q, - {Du)q^^^\ < ^ 

that is (|4.24p . The last case to consider is when k < je < h, that can be actually 
treated by a combination of the first two. It suffices to prove that the inequalities 
in display (|4.24p still hold. Indeed, the first inequality in (|4.24p follows exactly as 
in the second case. As for the second estimate in (|4.24p . let us remark that, as 
je G £, we have that \{Du)q.^ \ < Ae/50. On the other hand, we can use the first 
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case k < h < je with h — je, thereby obtaming \{Du)q^^ — (Du)q^| < Ae/50 and 
therefore the second inequaUty in (|4.24p follows via triangle inequality. 

5. Proof of Theorem 11.51 

The proof of Theorem 11.51 is a consequence of a few simple observations once 
Lemma [277l is at disposal and the sequence of comparison solutions Wj is introduced. 
Let's start by Theorem ll.il Going back to Section [2?2l after having introduced the 
maps {wj} in Lemma |2. 71 wc can introduce the maps {vj} exactly as in (j2.5|) . It is 
now easy to see that with this definition all the properties of the maps {wj} and 
{vj} described in Sections 12.2112.41 and used in the proof of Theorem 1 1 . 1 1 hold . and 
especially Lemma 12.4112.61 The only difference, which stems from the right hand 
sides in the inequalities in Lemma 12.71 is that the quantities iMKLQjJpar) appear 
instead of this is anyway irrelevant in the context of Theorem ll.il in view 

of (|3.13p and of the first equality in p.l5p . The proof now follows exactly as in the 
finite energy case. As a consequence, all the corollaries of Theorem ll.il starting by 
Theorem ll.21 hold for SOLA as well. Next, when passing to Theorems 11.31 and 11.41 
the proofs remain completely the same upon using Lemma 1 2. 71 note that this time 
the appearance of |M|([QjJpar) instead of \fJ-\{Qj) gives no problem at all since all 
the proofs arc based on the use of smallness conditions as (|4.5p and (|4.18|) . 

Remark 5.1. Definition [T] of SOLA is quite natural, and is motivated by the 
standard way of approximating measures with bounded functions in the weak-* 
convergence, via parabolic smoothing, that is, using mollifiers that, although acting 
uniformly in space, act backward in time. The main difference with the standard 
elliptic case is that, instead of getting that 

limsup \fih\iQ) < ImKQ) , 

h 

that is an inequality involving the full closure of Q, we get (|1.23p so that the upper 
part of a cylinder does not play any role in the approximation. This is the advantage 
that allows to pass to the limit easily in the pointwsc potential estimates. Moreover, 
this in accordance to the fact that, when dealing with evolutionary equations, the 
behavior at a certain instant of the solution only depends on what happened at past 
times, but not on what happens at that instant. Let us briefiy recall the procedure. 
One fixes a family of smooth mollifiers (approximation of the identity) {(f>h} with 
(j>h ■= h"-(j>{x/h) with (f> G C^{Bi), (f> > and ||0||li = 1- Similarly, we consider 
another family of smooth mollifiers {(l)h}, this time in one variable: 4'h ■= h(f>{x/h) 
with 4) e C(f ((-1, 1)) and ||^||ii = 1. We then define 



and notice that the mollification in the time variable is backward. In such a way we 
obtain a weakly* convergent sequence (in the sense of measures) and also (|1.23p . 
Such a sequence can be for instance used in [4] to derive the corresponding existence 
theorems. The point that we want to stress here is that, in order to give a more 
suitable definition of SOLA, it is very often not sufficient to take any approximation 
of the measure n via weak* convergence (something that is for instance sufficient 
when deriving global estimates). Instead, a more careful way of approximating 
measures, tailored to both the geometry of the problem in question and the degree 
of fine properties of solutions one wants to derive, must be adopted. 
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